EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-PH-EP /2007-026 

12 July 2007 

Measurement of the Mass and 
Width of the W Boson in e~^e~ 
Collisions at = 161 - 209 GeV 

DELPHI Collaboration 



Abstract 



A measurement of the W boson mass and width has been performed by the 
Delphi collaboration using the data collected during the full LEP2 programme 
(1996-2000). The data sample has an integrated luminosity of 660 pb~^ and 
was collected over a range of centre-of-mass energies from 161 to 209 GeV. 
Results are obtained by applying the method of direct reconstruction of the 
mass of the W from its decay products in both the W+W~ ^i^^qq' and 
W'^'W" —>■ qq'qq' channels. The W mass result for the combined data set is 

Mw = 80.336 ±0.055(Stat.)±0.028(Syst.)±0.025(FSl)±0.009(LEP) GeV/c^ 

where FSI represents the uncertainty due to final state interaction effects in 
the qq'qq' channel, and LEP represents that arising from the knowledge of the 
collision energy of the accelerator. The combined value for the W width is 

Fw = 2.404 ± 0.140(Stat.) ± 0.077(Syst.) ± 0.065(FSI) GeV/cl 

These results supersede all values previously published by the DELPHI collab- 
oration. 
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1 Introduction 

The measurement of the W boson mass can be used, in combination with other elec- 
troweak data, to test the vahdity of the Standard Model and obtain estimates of its 
fundamental parameters. In particular the measurement is sensitive, through loop cor- 
rections, to the masses of the top quark and the Higgs boson. 

The W boson mass and width results presented in this paper are obtained from data 
recorded by the Delphi experiment during the 1996-2000 operation of the Lep Collider, 
known as the Lep2 period. This corresponds to a total of 660 pb~^ collected over a range 
of centre-of-mass energies: ^/s — 161 — 209 GeV. 

Initially, data were recorded close to the W+W" pair production threshold. At this 
energy the W+W~ cross-section is sensitive to the W boson mass, Mw- Subsequently, 
Lep operated at higher centre-of-mass energies, where the e~^e~ — > W"'"W~ cross-section 
has little sensitivity to Mw For these data, which constitute the bulk of the Delphi data 
sample, Mw and the W boson width, T^^, are measured through the direct reconstruction 
of the W boson's invariant mass from the observed jets and leptons. The analysis is 
performed on the final states in which both W bosons in the event decay hadronically 
(W^W~ qq'qq' or fully-hadronic) and in which one W boson decays hadronically while 
the other decays leptonically (W+W" —>■ iP^qc^ or semi-leptonic). 

The Mw analyses of the relatively small quantity of data (~ 20 pb"^) collected during 
1996 at centre-of-mass energies of 161 and 172 GeV were published in [1,2]. These data 
are not reanalysed in this paper but are discussed in sections 7.1.1 and 7.1.2 and included 
in the final Mw combination. 

The data recorded during 1997 and 1998 at y/s = 183 and 189 GeV have also been 
the subject of previous Delphi pubhcations [3,4]. These data have been reprocessed 
and are reanalysed in this paper; the results given here supersede those in the previous 
publications. Results on the data collected during the final two years of Lep operation 
are published here for the first time. The data quality, simulation samples and analysis 
techniques have all been improved with respect to those used in previous Delphi publica- 
tions. The W mass and width have also been determined by the other Lep collaborations 
[5] and at hadron colliders [6]. 

The results on the W mass, Mw, and width, Fw, presented below correspond to a 
definition based on a Breit-Wigner denominator with an s-dependent width, | (s — Mw^) + 
isFw/MwI- 

After these introductory remarks, the paper starts in section 2 by describing the 
Lep accelerator and the determination of its collison energy. A brief description of the 
Delphi detector is provided as section 3. This is followed by section 4 which presents 
the properties of the data sample and of the Monte Carlo simulation samples used in the 

analysis. 

The analysis method is presented in section 5, first for W"'"W~ — > ^i^^qq' events, then 
for W+W" —>■ qq'qq' events. The text describes how the events are selected and the mass 
and width estimated from Mw- and Fw-dependent likelihood functions. The potential 
sources of systematic uncertainty are considered in section 6. These include: inaccuracies 
in the modelling of the detector; uncertainties on the background; uncertainties on the 
effects of radiative corrections; understanding of the hadronisation of the W boson jets; 
possible cross-talk between two hadronically decaying W bosons, the effects of which 
the qq'qq' Mw analysis has been specifically designed to minimise; and uncertainty on 
the Lep centre-of-mass energy determination. The paper concludes in section 7 with a 
presentation of the results and their combination. 
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2 LEP Characteristics 

2.1 Accelerator Operation 

The Lep2 programme began in 1996 when the colhsion energy of the beams was first 
ramped to the W+W" production threshold of 161 GeV and approximately 10 pb~^ of 
integrated luminosity was collected by each experiment. Later in that year Lep was run 
at 172 GeV and a dataset of similar size was accumulated. In each of the four subsequent 
years of operation the collision energy was raised to successively higher values, and the 
accelerator performance improved such that almost half the integrated luminosity was 
dehvered at nominal colhsion energies of 200 GeV and above. The main motivation for 
this programme was to improve the sensitivity of the search for the Higgs boson and 
other new particles. The step-by-stcp nature of the energy increase was dictated by the 
evolving capabilities of the radio frequency (RF) accelerating system. 

During normal operation the machine would be filled with 4 electron and 4 positron 
bunches at Ebeam ~ 22 GeV, and the beams then ramped to physics energy, at which 
point they would be steered into collision and experimental data taking begun. The 
fill would last until the beam currents fell below a useful level, or an RF cavity trip 
precipitated loss of beam. The mean fill lengths ranged from 5 hours in 1996 to 2 hours 
in 1999. After de-Gaussing the magnets the cycle would be repeated. 

In 2000, the operation was modified in order to optimise still further the high energy 
reach of Lep. Fills were started at a beam energy safely within the capabilities of the 
RF system. When the beam currents had decayed significantly, typically after an hour, 
the dipoles were ramped and luminosity delivered at a higher energy. This procedure was 
repeated until the energy was at the limit of the RF, and data taken until the beam was 
lost through a klystron trip. These mini-ramps lasted less than a minute, and varied in 
step size with a mean value of 600 MeV. The luminosity in 2000 therefore was delivered 
through a near-continuum of collision energies between 201 and 209 GeV. 

In addition to the high energy running, a number of fills each year were performed at 
the Z resonance. This was to provide calibration data for the experiments. Finally, several 
fills were devoted to energy calibration activities, most notably resonant depolarisation 
(RDP), spectrometer and Qg measurements (see below for further details). 

The machine optics which were used for physics operation and for RDP measurements 
evolved throughout the programme in order to optimise the luminosity at each energy 
point. Certain optics enhanced the build-up of polarisation, and thus were favoured for 
RDP measurements. The optics infiuence Ebeam in several ways, and are accounted for 
in the energy model, full details of which are available in [7] . 

2.2 The LEP Energy Model 

A precise measurement of the Lep beam energy, and thus the centre-of-mass energy, is 
a crucial ingredient in the determination of the W mass as it sets the overall energy scale. 
The absolute energy scale of Lep is set by the technique of RDP, which is accurate to 
better than 1 MeV. This technique allowed very precise measurements of the mass and 
width of the Z boson to be made at LepI. However, this technique is only possible for 
beam energies between about 41 and 61 GeV. The Lep2 energy scale is set mainly by the 
nuclear magnetic resonance (NMR) model. This makes use of 16 NMR probes, positioned 
in selected dipoles, which were used to obtain local measurements of the bending field. 
These probes thus sample the total bending field, which is the primary component in 
determining the beam energy. Onto this must be added time-dependent corrections 
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coming from other sources. These include effects from earth tides, beam orbit corrections, 
changes in the RF frequency, and other smaller effects. Details of all these can be found 
in [7] . Using this Lep Energy Model, the Lep Energy group provided Delphi with an 
estimate of the centre-of-mass energy at the start of each fill and thereafter in intervals 
of 15 minutes. For the year 2000 the values before and after the mini-ramps were also 
supplied. No data are used which are taken during the mini-ramps, as the energy is not 
accurately known during these periods. 

The main assumption which is made in the Lep Energy Model is that the beam energy 
scales linearly with the readings of the NMR probes. This assumption of linearity has 
been tested by three different methods: 

1) Flux Loop. Each dipolc magnet of Lep is equipped with a single-turn flux loop. 
Measurements are made for a series of dipole magnet currents, which correspond 
roughly to the operating beam energies of Lep2. This allows the change in flux 
over almost the entire LEP dipole field to be measured as the machine is ramped in 
dedicated experiments. This change in flux can be compared with the local bending 
flcld measurements of the NMR probes. The Flux Loop is calibrated against the 
Lep energy model in the range 41-61 GeV, using the NMR coefficients determined 
from RDP. The measurements from the Flux Loop in the high energy regime (up 
to 106 GeV beam energy) are then compared to those from the Lep Energy Model. 
The Flux Loop measurements were made in all years of Lep2 running. 

2) Spectrometer Magnet. In 1999 a special steel Spectrometer Magnet, equipped with 
three beam position monitors to measure the beam position both on entry and 
exit from the magnet, was installed in the Lep ring. The magnetic field of this 
magnet was carefully mapped before and after installation in the Lep ring. All these 
measurements were very compatible. The beam energy is determined by measuring 
the bending angle of the beam in passing through the dipole magnet. The device 
was calibrated against RDP in the 41-61 GcV region and the Spectrometer results 
were compared to the Lep Energy Model at beam energies of 70 and 92 GeV. 

3) Qs versus Vrf- The synchrotron tune Qs can be expressed as a function of the 
beam energy and the total RF voltage, Vrf, plus some additional small corrections. 
By measuring Qs as a function of the total RF voltage the beam energy can be 
determined. These measurements were performed in 1998-2000, at beam energies 
from 80 to 91 GeV. Again the measurements were normalised against RDP in the 
region 41-61 GeV, and compared to the Lep Energy Model at Lep2 energies. 

The three methods are in good agreement, both with each other and the Lep Energy 
Model. Based on these comparisons a small energy offset compared to the Lep Energy 
Model was supplied for each of the 10 beam energies used in Lep2. This offset is always 
smaller than 2 MeV. The estimated centrc-of-mass energy uncertainties range between 
20 and 40 McV and are discussed further in section 6.8 . 

The Lep centre-of-mass energy has also been determined by the Lep collaborations 
using Lep2 events containing on-shell Z bosons and photons (radiative return to the Z 
events) [8,9]. The Delphi analysis measured the average difference between the centre- 
of-mass energy from radiative return events in the e~ fi'^ fi~ {'j) and e'^e~— ^55(7) 
channels and the energy reported by the Lep Energy working group. 



= +0.073 ± 0.094:{Stat.) ± 0.065{Syst.) GeV. 
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Thus the Delphi result, relying on similar reconstruction procedures to those de- 
scribed in this paper, is in agreement with the values reported by the Lep Energy working 
group. 

3 Detector Description 

The Delphi detector [10] was upgraded for Lep2. Changes were made to the sub- 
detectors, the trigger system, the run control and the algorithms used in the offline 
reconstruction of tracks, which improved the performance compared to the earher LepI 
period. 

The major change was the inclusion of the Very Forward Tracker (VFT) [11], which 
extended the coverage of the innermost silicon tracker out to 11° < < 169°^. Together 
with improved tracking algorithms and alignment and calibration procedures optimised 
for Lep2, these changes led to an improved track reconstruction efficiency in the forward 
regions of Delphi. 

Changes were made to the electronics of the trigger and timing system which improved 
the stability of the running during data taking. The trigger conditions were optimised for 
Lep2 running, to give high efficiency for Standard Model two- and four-fermion processes 
and also to give sensitivity for events which may be signatures of new physics. In addition, 
improvements were made to the operation of the detector during the Lep cycle, to prepare 
the detector for data taking at the very start of stable collisions of the e^e" beams, and 
to respond to adverse background from Lep were they to arise. These changes led to 
an overall improvement of ~ 10% in the efficiency for collecting the delivered luminosity 
from ~ 85% in 1995, before the start of Lep2, to ~ 95% at the end in 2000. 

During the operation of the DELPHI detector in 2000 one of the 12 sectors of the 
central tracking chamber, the TPC, failed. After the P*^ September 2000 it was not 
possible to detect the tracks left by charged particles inside the broken sector. The data 
affected correspond to ~ 1/4 of the total dataset of the year 2000. Nevertheless, the 
redundancy of the tracking system of DELPHI meant that tracks passing through the 
sector could still be reconstructed from signals in any of the other tracking detectors. 
A modified track reconstruction algorithm was used in this sector, which included space 
points reconstructed in the Barrel RICH detector. As a result, the track reconstruction 
efficiency was only slightly reduced in the region covered by the broken sector, but the 
track parameter resolutions were degraded compared with the data taken prior to the 
failure of this sector. 

4 Data and Simulation Samples 
4.1 Data 

The W mass and width are measured in this paper with the data samples collected 
during the 1996-2000 operation of the Lep Collider. A summary of the available data 
samples is reported in table 1, where the luminosity- weighted centre-of-mass energies and 
the amount of data collected at each energy are shown. The luminosity is determined 
from Bhabha scattering measurements making use of the very forward electromagnetic 
calorimetry [12]. The total integrated luminosity for the Lep2 period corresponds to 

^The Delphi coordinate system is right-handed with the z-axis coUinear with the incoming electron beam, and the x 
axis pointing to the centre of the Lep accelerator. The radius in the xy plane is denoted R and 6 is used to represent the 
polar angle to the z axis. 
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approximately 660 pb~^ The integrated luminosities used for the different selections 
correspond to those data for which all elements of the detector essential to each specific 
analysis were fully functional. The additional requirements on, for example, the status 
of the calorimetry and the muon chambers mean that the integrated luminosity of the 
semi-leptonic analysis is slightly less that that of the hadronic datasct. 

All the data taken from the year 1997 onwards have been reprocessed with an improved 
reconstruction code, and the analyses on these data are updated with respect to the 
previously published ones and supersede them. The data taken in 1996 have not been 
reanalysed; the results from this year are taken from the previous publications with minor 
revisions as reported in section 7. 

In addition to these data taken above the W+W^-pair production threshold, data were 
also recorded during this period at the Z peak. These samples, containing a total of over 
0.5 million collected Z decays, were taken each year typically at the start and end of the 
data taking periods. These Z peak samples were used extensively in the alignment and 
calibration of the detector and are used in many of the systematic uncertainty studies 
reported in section 6. 



Year 


/^-AA'oighted \/7i (GeV) 


Htulroiiic Int. C (pb 


Loptonic Int. C (pb ^) 


1996 


161.31 


10.1 


10.1 




172.14 


10.1 


10.1 


1997 


182.65 


52.5 


51.8 


1998 


188.63 


154.4 


152.5 


1999 


191.58 


25.2 


24.4 




195.51 


76.1 


74.6 




199.51 


82.8 


81.6 




201.64 


40.3 


40.2 


2000 


205.86 


218.4 


215.9 



Table 1: Luminosity- weighted centre-of-mass energies and integrated luminosities in 
the Lep2 data taking period. The hadronic integrated luminosity is used for the ful- 
ly-hadronic channel, the leptonic one is used for the semi-leptonic channels. 



4.2 Simulation 

The response of the detector to various physical processes was described using the 
simulation program DELSIM [10], which includes modelling of the resolution, granular- 
ity and efficiency of the detector components. In addition, detector correction factors, 
described in section 6, were included to improve the description of jets, electrons and 
muons. To allow use of the data taken after the 1®* September in 2000, samples of events 
were simulated dropping information from the broken sector of the TPC. A variety of 
event generators were used to describe all the physics processes relevant for the analy- 
sis. W'^W" events and all other four-fermion processes were simulated with the program 
described in [13], based on the WPHACT 2.0 generator [14] interfaced with PYTHIA 6.156 
[15] to describe quark hadronisation and TAUOLA 2.6 [16] to model r leptons decays. The 
most recent 0{a) electroweak radiative corrections in the so-called Double Pole Approx- 
imation (DPA) were included in the generation of the signal via weights computed by 
YFSWW 3.1.16 [17], and the treatment of initial state radiation (ISR) of this calculation was 
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adopted. The photon radiation from final state leptons was computed with PHOTOS 2.5 
[18]. For systematic studies the alternative hadronisation descriptions implemented in 
ARIADNE 4.08 [19] and HERWIG 6.2 [20] were also used. All the hadronisation models were 
tuned on the Delphi Z peak data [21]. 

The background process e^e~ QQil) was simulated with KK 4.14 [22] interfaced 
with PYTHIA 6.156 for the hadronisation description. The two-photon events giving rise 
to those e'^e~qq final states not described in the four-fermion generation above were pro- 
duced with PYTHIA 6.143 as discussed in [13]. The contribution from all other background 
processes was negligible. 

The simulated integrated luminosity used for the analysis was about a factor 350 
higher than for the real data collected for 4-fermion processes, about a factor 60 higher 
for 2-fermion final states and about 3.5 times greater for e'^e''qq two-photon final states 
(those not already included in the 4-fermion simulation). 

5 Analysis Method 

The measurement of Mw and of Fw are performed on samples of W"'"W~ — > iUeqc^ 
and W+W" qq'qq' events; these two channels are discussed in turn below. The 
reconstruction of events where both Ws decay leptonically has very limited sensitivity to 
the W mass and width, as they contain at least two undetected neutrinos, and hence are 
not used in this analysis. 

The first stage in the analysis is to select events from these decay channels, using 
either a neural network or a sequential cut-based approach. In some channels, after 
preliminary cuts, the probability is assessed for each event of how W'*"W~-like it is and 
a corresponding weight is applied in the analysis. 

The resolution of the kinematic information extracted from the observed particles in 
the event can be improved by applying energy and momentum conservation constraints to 
the event; this is discussed in section 5.1. In the fuUy-hadronic channel the jet directions 
used as the input to the kinematic fit are also assessed excluding particles from the inter- 
jet regions. This alternative approach reduces the sensitivity of the W mass analysis to 
final state interaction systematics and is discussed in section 5.3.2. 

The next stage in the analysis is to produce a likelihood function expressing the relative 
probability of observing an event as a function of Mw and Fw- The likelihood functions 
used below depend not only on the reconstructed W mass of the event but make use 
of other event characteristics to assess the relative weight and resolution of each event. 
These likehhood functions are then cahbrated against simulated events. 

The W mass and width are then extracted by maximising the combined likelihood 
function of the full observed dataset. 

5.1 Application of Kinematic Constraints to Event Reconstruc- 
tion 

The event-by-event uncertainty on the centre-of-mass energy, i. e. the energy spread, at 
Lep is typically 0.1%, while the overall momentum and energy resolution of the observed 
final state is about 10%. Hence, the precise knowledge of the kinematics in the initial state 
can be used to significantly improve the reconstructed kinematic information obtained 
from the clustered jets and observed leptons in the final state. This is accomphshed by 



7 



means of a fit based on the four constraints from the conservation laws of energy and 
momentum. 

The reconstructed jets and leptons of the event may be associated with one of the two 
hypothesised W bosons in the event. A fifth constraint may then be apphed to the event 
by assigning equal masses to these W boson candidates. As the decay width of the 
bosons is finite, this constraint is non-physical. However, as the event mass resolution 
and 2 GeV/ W width are of comparable magnitude in practice this constraint provides 
a useful approximation. It is of particular use in the semi-leptonic decay channels where, 
after applying the four-constraints, the event mass resolution is still larger than the W 
width and, due to the unseen neutrino, the two fitted masses are strongly anticorrelated. 
However, in the fuUy-hadronic decay channel the mass resolution after the four-constraint 
fit is better and the correlation is less; hence more information is available in the two 
four-constraint masses than the combined five-constraint event mass. 

Pcirameterisation of Jets and Leptons 

Each fitted object, jet or lepton, is described by three parameters. Muons are described 

by their measured momenta and their polar and azimuthal angles. The uncertainties on 
these parameters arc obtained directly from the track fit. Electrons arc characterized 
by their measured energies and their detected angular position in the electromagnetic 
calorimeters. The energy uncertainties are obtained from parameterisations of the re- 
sponses of the electromagnetic calorimeters, which were tuned to the responses found 
in Bhabha and Compton scattering events. The angular uncertainties were determined 
from the detector granularity and were significant only for the forward electromagnetic 
calorimeter. In semi-leptonic events, the neutrino momentum vector is considered as un- 
known, which leads to a reduction by three in the number of effective constraints in the 
kinematic fit. 




Figure 1: Parameterisation used for jets in the constrained fit, as explained in the text 
and equation 1. 
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Each fitted jet momentum is projected onto a set of axes with one component 
parallel to the measured jet momentum p^™" and two transverse components, "p^ and p^*^, 

each normalized in magnitude to 1 GeV/ c. In this coordinate system can be described 
by three parameters a^, hj and cf 

p/ = e'^^p,.'" + 6,p/ + c,p/, (1) 

where each component is shown in figure 1. The measured jet energy E.^ is rescaled 
with the same factor e"^ as the jet momentum. The exponential parameterisation e"^ of 
the factor in front of p^"* makes the fit more stable and results in uncertainties which 
have a more Gaussian distribution. The values of the parameters are determined by 
performing a constrained fit, while the transverse directions are given by the eigenvectors 
of the momentum tensor described below. 

Form of 

The algorithm minimizes a x^, defined for fuUy-hadronic events as: 

■jets I \2 L 2 2 

j=l % % ^c, 

while forcing the fitted event to obey the constraints. The appropriate terms are included 
in the x^ for events with a leptonic W decay. The expected energy loss parameter ao 
and the energy spread parameter a^, together with the parameters ai,. and ctc^, are 
parameterised as functions of the jet polar angles. 

Jet Error Pcirameterisation 

The jet error parameters, ao, (7^ , (Jbj and were obtained from a study of hadronic Z 
events. Hadronic Z events with a two-jet topology were selected from the Z calibration run 
data or from the corresponding Monte Carlo simulation. The reconstructed jet energies 
were compared with the beam energy. In general an energy loss of around 10% was 
observed for jets in the barrel region of the detector while this increased to 15% in the 
forward regions. A good agreement between the data and simulation was found. The 
energy loss increases if the event jet topology becomes less two-jet like, resulting in energy 
losses of around 15% for the barrel region and up to 35% in the forwards regions. 

The uncertainties on the jet parameters for the first stage of the fit were determined 
from this study as a function of the polar angle of the jet. However, a dependence of 
these parameters on the properties of the individual jets has also been observed. 

Jet Breadth 

The dependence of the uncertainties on the individual jet properties is included in a 

second stage of the fit, where the parameterisation of the transverse momentum uncer- 
tainties depends upon the breadth of the jet. This breadth is calculated by projecting 
the momenta of all particles in the jet on to the plane transverse to the jet axis. From 
these projections a two dimensional momentum tensor is created: 

= (3) 

k 

where and are the two components of the projection of the momentum of particle 
k in the transverse plane. The normalized eigenvectors of the tensor, p-^ and p-^, reflect 
the directions where the jet is broadest and slimmest. The corresponding eigenvalues are 
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Bi, and B^. By comparing the resulting jet energies from the first stage of the fit with 
the measured ones, an estimate is made of how much energy remained undetected in the 
jet, referred to as Ej^rniss- The uncertainties on the jet breadths were then parametrised 
as a function of the eigenvalues, the measured jet energy and the missing energy Ej^rniss- 

Use of 

The of the resulting fit is a function of the collection of jet parameters (a^, hj, cj) 
and lepton parameters. The jets and leptons are paired appropriately to each W boson 
decay and constraints applied. The total is then minimized by an iterative procedure 
using Lagrange multipliers for the constraints. 

Events for which the of the fit is larger than the number of degrees of freedom for the 
fit, NDF, had their errors scaled by a factor of ^Jx^/NDF in order to take non-Gaussian 
resolution effects into account. 

In the semi-lcptonic analysis described in section 5.2.3 the value of the best fit mass 
from the x^ minimum and the error on this mass is used for each event. In the fuUy- 
hadronic analysis described in section 5.3.3 each event uses the x^ distribution as a 
function of the masses of the two W bosons in the event. 

5.2 Semi-Leptonic Decay Channel 

The W+W- iUeqc^ events constitute 44% of all W+W" decays. The W+W" event 

candidates are classified according to their leptons and their selection is performed using 
a neural network. An event W mass is reconstructed in a kinematic fit, by imposing 
momentum conservation, the measured centre-of-mass energy and equality of the leptonic 
and hadronic decay W masses. An estimate of the mass resolution in each individual event 
is also obtained from the kinematic fit and an estimate of the event purity is obtained 
from the neural network output; these quantities are both used in producing the likelihood 
function from which Mw and Fw are determined. 

5.2.1 Event Selection 

Events are selected from the recorded data sample requiring that all detectors essential 
for this measurement were fully efficient: these comprise the central tracking detectors 
and the electromagnetic calorimeters. The data recorded during the period with a dam- 
aged sector of the TPC are also used with matching simulation samples produced. The 
corresponding integrated luminosities, at each centre-of-mass energy, are given in table 1. 

Events containing at least three charged particle tracks and with a visible mass greater 
than 20 GeV/ are considered for analysis. Events containing lepton candidates are then 
identified in this sample, either by direct lepton identification (electrons and muons), or by 
clustering the events into a three-jet configuration and selecting the jet with the lowest 
charged multiplicity as the tau candidate. At this stage, events can be considered as 
candidates in multiple channels. 

Electron and Muon Identification 

Charged particles are identified as muons if they are associated with a hit in the muon 
chambers, or have an energy deposit in the hadron calorimeter that is consistent with 
a minimum ionising particle. Muon identification is performed in the polar angle range 
between 10° and 170°. Muons with an unambiguous association [10] with the hits in 
the muon chambers, or with a loose association in addition to a good pattern in the 
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hadron calorimeter are classified as good candidates, with the remainder being classified 
as possible candidates. 

Electron identification is performed in the polar angle range between 15° and 165° by 
selecting charged particles with a characteristic energy deposition in the electromagnetic 
calorimeters. In the central region of the detector, covered by the HPC electromagnetic 
calorimeter, the electron selection followed the criteria described in [10] for candidates 
below 30 GeV. This selection is based on a neural network using the electron energy 
to momentum ratio (E/p), the spatial matching between the extrapolated track and 
the shower, the shower shape and the track energy loss per unit path length in the 
TPC (dE/dx) as the discriminating variables. Above 30 GeV, a simplified selection 
is adopted, the main deposit associated with a charged particle track is identified and 
the surrounding electromagnetic showers are clustered into this electron candidate. Only 
candidates with E/p greater than 0.5 are used. In the polar angle region corresponding to 
the forward electromagnetic calorimeter acceptance, below 36° and above 144°, electron 
candidates are selected from among the calorimetric shower clusters. Only clusters with 
an energy above 8 GeV and which could be geometrically associated to extrapolated 
charged particle tracks arc used. The electron candidates arc separated into categories 
of good and possible candidates based on the quality of the track associated with the 
electron. The association of vertex detector hits to the track is a primary criterion used 
in assessing the track quality. 

Tau reconstruction 

As mentioned above, tau candidate events arc clustered into a three-jet configuration 
using the LUCLUS [23] algorithm. Tracks at large angle (more than 40° from the nearest 
jet axis) or which contribute a large mass to the jet they belong to (AM bigger than 
3.5 GeV/c^) are removed from the tau candidate. As the tau lepton predominantly 
decays into a final state with one or three charged particles, with few neutrals, a pseudo- 
multiplicity defined as the sum of the charged multiplicity and one quarter of the neutral 
multiplicity is used and the jet with the lowest pseudo-multiplicity is chosen as the tau 
candidate. Then a further cleaning is applied on this tau candidate : tracks at more than 
20° from the tau axis, or which contribute a large mass (AM bigger than 2.5 GeV/ c^) are 
removed from the tau candidate. Only tau candidates containing between one and four 
charged particle tracks after this cleaning, and with a polar angle between 15° and 165° 
are kept. Two classes of events are then defined, those with only one charged particle 
track, and all others. 

Event Reconstruction and Pre-selection 

After the lepton identification is performed, the events are reconstructed as the lepton 
and a two or three jet system. Pre-selection cuts are then applied. 

All tracks not associated to the lepton are clustered using the LUCLUS algorithm. These 
jet tracks in semi-leptonic electron and muon decay channel events are clustered with 
djoin = 7.5 GeV/ c,where djoin is a measure of the clusterisation scale used inside LUCLUS. 
If more than three jets are obtained the tracks are forced into a three-jet configuration. 
This procedure correctly treats events with hard gluon radiation (the proportion of three- 
jet events is about 20%). In semi-leptonic tau decay events the tracks not associated to 
the tau candidate are forced into a two-jet configuration. 

A set of pre-selection cuts is then applied. First, a common set of criteria is applied 
to the system of jets: 

• visible mass greater than 30 GeV/c^; 
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• at least five charged particle tracks, with at least two with momentum transverse to 
the beam greater than 1.5 GeV/c and compatible with the primary vertex (impact 
parameter in it! < 0.15 cm and in z < 0.4 cm); 

• no electromagnetic cluster with an energy bigger than 50 GeV. 

Then, for electron and muon semi-leptonic decay channel events, the following addi- 
tional cuts are used: 

• energy of the lepton bigger than 20 GeV; 

• if there is another isolated lepton of the same flavour and opposite charge, the 
event acollinearity should be bigger than 25°. The acollinearity used here is that 
between the two 'jets' when forcing the event into a two-jet (including the lepton) 
configuration. 

Further cuts are made for electron decay channel events: 

• missing transverse momentum should be greater than 8 GeV/c; 

• the cut on missing transverse momentum is increased to 12 GeV/c for electron 
candidates in the 'possible' class; 

• angle between the lepton and the nearest jet greater than 15°. 

The cuts specific to the muon decay channel events are: 

• angle between the lepton and the nearest jet greater than 15° in the case of 'possible' 
class muons; 

• angle between the missing momentum and the beam axis greater than 10° for muon 

candidates in the 'possible' class. 

While for tau decay channel events, the cuts applied are: 

• visible hadronic mass smaller than 130 GeV/ c^; 

• energy of the tau greater than 5 GeV; 

• fraction of energy of the tau associated to charged tracks greater than 5%; 

• at least one of the charged particle tracks from the tau must have a vertex detector 
hit; 

• angle between the tau and the nearest jet greater than 15°; 

• angle between the tau and the nearest charged particle greater than 10°; 

• missing transverse momentum greater than 8 GeV/ c; 

• the cut on missing transverse momentum is increased to 12 GeV/ c in the case of 
tau candidates with several charged particles. 

The semi-leptonic electron and muon events are then reconstructed using a constrained 
fit imposing conservation of four-momentum and equality of the two W masses in the 
event. As the energy of the tau lepton is unknown, due to the emission of at least one 
neutrino in its decay, the mass in the rz7,-qq' channel is entirely determined by the jet 
system and no improvement can be made from applying a constrained fit. 

Selection 

The event selection is based upon a multi- layer perceptron neural network [24]. The 
network has been optimised separately for the six classes of events (good and possible 
eI7eqq', good and possible /xl/^qq', and Ti^T-qq' candidates containing either only one or 
several charged particles). 

The choice of the variables used in the neural networks is a compromise between 
their independence from the W mass and their discriminant power. The number of 
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Figure 2: The output of the neural network used for the selection of the semi-leptonic 
channels for the data sample recorded at y/s = 183 — 209 GeV. The data are indicated 
by the data points with error bars. The histograms show the signal and background 
simulation contributions normalised to the integrated luminosity of the data sample. 
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input-hidden-output nodes were 12-8-1, 11-7-1 and 17-12-1 for the e, and r channels 
respectively. The detailed list of variables is given below. The network has been tuned 
on samples of signal and background simulation events, and examples of the distribution 
of the neural network output value are shown in figure 2. The applied selection cut is at 
0.40, 0.50 and 0.35 for the e, ^ and r channels respectively, independent of the centre- 
of-mass energy. Any discrepancy in the background rate between data and simulation is 
accounted for in the systematic uncertainty applied. 

The event selection procedure ensures that the events are only selected in one of the 
channels: events that pass the chosen cut in the muon channel are selected, the remaining 
events are considered as electron channel candidates and, if they are again rejected, are 
then analysed under the tau channel hypothesis. This ordering follows the hierarchy of 
purities in these channels (and is not dependent on the good or possible lepton classes). 
After applying the cut on the network output the selection performance is as shown in 
table 2. As an example, the global efficiencies for CC03 events arc 79.8%, 89.8% and 59.3% 
respectively for the eI7eqq', nV^qq^ and rl^T-qq' events in the data taken at y/s = 189 GeV. 
These numbers are integrated over all event selections as there is a non-negligible cross- 
contamination of events in the event selections {e.g. ez7eqq' event selected by the rz/T-qq' 
selection) which still add useful information in the W mass and width fits. Here CC03 
refers to the three charged current processes producing the W"'"W^ state for which this 
analysis is intended: s-channel photon or Z production and t-channel exchange. 
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Data 


el/eqq' 


257.5 


10.5 


0.7 


9.3 


6.5 


284.5 


256 


/iz7^qq' 
rI7^qq' 


321.2 


10.2 


0.4 


1.1 


2.2 


335.1 
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198.2 


56.6 


3.5 


18.6 


10.9 


287.9 


294 


qq'qq' 




34.0 


1029.9 


341.6 


50.8 


1456.3 
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2000, 206 GeV 
















ez7cqq' 


373.9 


16.9 


1.0 


13.6 


11.4 


416.8 


395 


/iz7^qq' 


457.0 


14.8 


0.6 


1.7 


4.1 


478.2 


467 


Tl^T-qq' 


290.2 


87.6 


5.7 


22.3 


21.4 


427.2 


426 


qq'qq' 




40.6 


1514.5 


460.9 


107.8 


2123.8 


2134 


1997-2000 
















183-206 GeV 
















ez/cqq' 


1091.5 


47.7 


2.9 


39.9 


30.7 


1212.7 


1182 


^F^qq' 
r77^qq' 


1356.7 


43.3 


1.7 


15.2 


11.0 


1417.8 


1402 


849.3 


248.6 


16.0 


72.2 


55.6 


1241.6 


1270 


qq'qq' 




131.6 


4421.0 


1399.5 


269.8 


6222.0 


6446 



Table 2: Number of selected events in the decay channel event selections from the 1998 
and 2000 data samples and the combined 1997-2000 data sample, and the corresponding 
number of expected events from the simulation. The table is split into rows giving the 
results of each of the event selection routines. The primary-/ and other-/ z/^qq' columns 
relate to the nature of the semi-leptonic event selections e.g. for the el^eqq' selection the 
results are for the el^eqq' and (//i^^nqq' + rl/^qq') channels respectively. 

For each of the six classes of events, the fraction of semi-leptonic W+W" events in the 
sample has been extracted from simulation as a function of the neural network output: 
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this is referred to below as the event purity Pg- This feature is particularly useful for the 
tau selection, where the proportion of background events is highest. 

5.2.2 Veiriables used in the Selection Neural Networks 
Common Variables for all Leptonic Channels 

• Polar angle of the leptonic W (after applying the constrained fit); 

• angle of the charged lepton with respect to the direction of the leptonic W (in the 
W rest frame, and after the constrained fit); 

• polar angle of the lepton; 

• polar angle of the missing momentum vector; 

• angle between the lepton and the nearest jet; 

• angle between the lepton and the nearest charged hadron track (of energy greater 

than 1 GeV); 

• missing transverse momentum; 

• the invariant mass of the measured system of particles [25] - this is measured 
using planar kinematics, by forcing the event into 2 jets (using all particles in the 
event including the lepton) and assuming a photon is emitted down the beam pipe; 

• aplanarity (cosine of the angle between the lepton and the normal to the plane 
formed by the jets^); 

• acoUinearity (complement of the angle between the two "jets" when forcing the event 
into a two-jet configuration); 

• the minimum djom distance in the LUCLUS jet clusterisation algorithm between two 
jets in the final configuration, where the whole event (hadronic and leptonic system) 
is forced into three jets. This is known as djsaii. 

Additional Variable for the Electron Channel Only 

• Angle between the missing momentum and the nearest jet. 
Additional Vctriables for the Tau Channel Only 

• Angle between the missing momentum and the nearest jet; 

• fraction of the tau energy coming from charged particle tracks; 

• missing energy; 

• reconstructed tau energy; 

• reconstructed tau mass; 

• dj4aii, as dj3aii (see above) but with the final event configuration forced into four jets. 

5.2.3 Likelihood Function 

A likelihood function, £e(Mw,rw)i is evaluated for each selected event with a recon- 
structed mass in a defined range. The range was 67-91 GeV/ for the data collected in 
1997, 67-93 GeV/c^ for 1998, 67-95 GeV/c^ for 1999, and 67-97 GeV/c^ for 2000. The 
increase in range with rising centre-of-mass energy is to account for the increasing ISR 
tail. The hkelihood function is defined as follows: 

£e(Mw, Tw) = Pe ■ S"{m^'\ <T^'\ Mw, Tw) + (1 - Pe) ■ B{mf''), 

where Pg is the event purity, discussed above, S" is the signal function that describes the 
reconstructed mass distribution of the semi-leptonic W decays, and B is used to describe 

*^for thrcc-jcts events in the electron and muon channels, the jets-plane is the plane formed by the most energetic jet 
and the sum of the two others. 
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background processes. The reconstructed event mass m'^** and its estimated error a^^^ 
are both obtained from the constrained fit. The distribution of background events is 
extracted from simulation as a function of m^''^. 

The signal function S" is defined in terms of S and S' as discussed below. The function 
S relies on the convolution of three components, using x and m as the dummy integration 
variables: 

5(m^**,cr^'*|Mw,rw) = (4) 
/(f^^*" dm G[mf'^ - m, a^'^] /q dx PS{m{l - x)) BW[m{l - a;)|Mw, Tw] Risr{x)- 
BW is a relativistic Breit-Wigner distribution representing the W mass distribution, 

Siy(m|Mw,rw) = To, (5) 

TT Mw (^2 _ + (m2^) 

and PS is a phase-space correction factor 



PS{m) 



Arm? 



! s 

The convolution with the Gaussian function G describes the detector resolution. The 
width of the Gaussian depends upon the reconstructed mass error obtained in the con- 
strained fit for that event. 

The ISR spectrum is parameterised as 

where x^ is the ratio of the photon energy to the ccntrc-of-mass energy and j3 is calculated 
from the electromagnetic coupling constant (a), the centre-of-mass energy squared (s) 
and the electron mass (mg): 

/? = — [ln(s/m^)-l]. 

TT 

Due to the constrained fit, a W produced at mass m will be reconstructed to a good 
approximation as m/(l — x^,) in the presence of an undetected ISR photon, giving a tail 
at high mass in the measured spectrum. This tail is well described by the integration on 
the photon spectrum in equation 4. 

The event selection contains a significant fraction of rl^T-qq' events in the electron and 
muon channel samples, and of el/cqq' events in the tau sample (see table 2). In the tau 
channel the mass of the event is determined from the jet system. The behaviour of true 
TUrOSl and eFeqq' events in this fit are found to be similar, and S" — S va. this channel. 
However, in the electron and muon channel samples the behaviour of the rl^T-qq' events is 
somewhat different to that of the eI7eqq', /iX^^qq' events. The rz/T-qq' events have a worse 
mass resolution and introduce a small negative offset on the mass. The fraction of tau 
events, which have been wrongly classified and are contained in the electron and muon 
channel samples, has been parameterised in bins of the lepton energy and the measured 
missing mass. This fraction Pt^ is then taken into account in the likelihood function for 
the electron and muon samples, by defining the signal function S" as 

S" ^{l-PTe)-S + PT,-S\ 

where <S" is analogous to S, but with the width of the Gaussian resolution function 
increased according to a factor determined from simulation studies. All remaining biases 
in the analysis due to using this approximate likelihood description are corrected for in 
the cahbration procedure as described in section 5.4. 
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5.3 Fully-Hadronic Decay Channel 

The W+W~ qq'qq' events constitute 46% of all W^W~ decays. The event masses 
can be reconstructed from the observed set of jets. The kinematics of the jets can be 
significantly over-constrained in a kinematic fit, improving the event mass resolution, by 
imposing momentum conservation and the measured centre-of-mass energy. The influence 
of the many ambiguities in the event reconstruction, which dilute the statistical informa- 
tion, is minimised by optimally weighting the different hypotheses in the likelihood fit of 
Mw or Fw- 

The dominant systematic error is due to the possible influence of flnal state interference 
effects between particles from the two decaying Ws. Reconstructing the jet directions 
using only the particles from the core of the jet reduces the possible effects of these final 
state interference effects. This technique and the mass estimator based on all observed 
particles are both discussed in section 5.3.2. 

5.3.1 Event Selection 

As in the semi-leptonic analysis, appropriate criteria were imposed on the functionality 
of the detector when selecting the data sample for analysis. The corresponding integrated 
luminosities, at each centre- of- mass energy, are given in table 1. 

The event selection can be separated into three stages. First a prc-selection is per- 
formed to reduce the data sample to events with a high multiplicity and high visible 
energy. In the second stage events with a four or five jet topology are retained. The 
observablcs on which the selection is made are chosen to be, to a good approximation, 
independent of the ccntrc-of-mass energy ^/s: the same selection criteria are used for 
all energies for the pre-selection and jet topology selection. The final stage of the event 
selection is to use the inter-jet angles and jet momenta to estimate the probability that 
this was a W"'"W^ — ^> qq'qq' event. 

The pre-selection cuts applied are: 

• the charged particle multiplicity should be larger than 13; 

• the total visible energy of the event must exceed 1.15 

• the scaled effective centre-of-mass energy ^ [25] is required to be equal to or larger 
than 0.8; 

• rejection of events tagged as likely to be containing b quarks [26]. 

The last criterion removes 7% of the remaining X^qqi^) and 18% of the remaining ZZ 
events, while changing the signal selection efficiency by less than 1%. The distributions 
of data and simulation events for the scaled effective centre-of-mass energy and combined 
b-tag variable are shown in figure 3; the cut on the combined b-tag variable retains all 
events below 2. 

The remaining events are then clustered using the DURHAM [27] jet clustering algorithm 
with a fixed ycut of 0.002. The jets obtained are required to have an invariant mass of 
greater than 1 GeV/ and contain at least three particles. If the jets do not meet these 
criteria or more than five jets are obtained, the clustering is continued to higher values 
of Ucut- Events which cannot be clustered into either four or five jets that fulfill these 
criteria are rejected. The initial ycut value of this procedure was optimised for maximal 
sensitivity to Mw and results in a sample of approximately 50% four and 50% five jet 
events. 

The jets obtained from this procedure are then used in a constrained fit, described in 
section 5.1, where momentum conservation and the measured centre-of-mass energy are 



17 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 -4 -2 2 4 6 8 10 12 

VsVVs Combined b-tag variable 



Figure 3: The distribution of two event selection variables for candidate qq'qq' events from 
the full Lep2 data sample and the corresponding simulation samples. The left hand plot 
shows the scaled effective centre-of-mass energy, the right hand plot the combined b-tag 
variable. The distributions are shown after the cuts on all other pre-selection variables 
have been applied. 

enforced. From the fitted jets a topological observable, Dpur, was formed to discriminate 
between signal events and Z—>qq events with hard gluon radiation: 

where and are the smallest and second smallest fitted jet energies and 9{j^ and 
6{j^ are the smallest and second smallest fitted inter-jet angles. The expected fraction of 
qq'qq' events (W^W~ or ZZ) in the selected sample, the event purity P^-^ , is parame- 
terised as a function of this variable. This fraction of qq'qq' events, i.e. doubly-resonant 
events rather than just W"'"W~ events, is used in the theoretical distribution function de- 
scribed below. Events with an estimated purity below 25% are rejected. The distribution 
of the Dpur observable is shown in figure 4 for both the 4 and 5 jet topology events, and 
the numbers of selected events are given in table 2. An excess of data events over the 
expected number of simulation events was observed. 

5.3.2 Cone Jet Reconstruction 

The largest contribution to the systematic uncertainty in the fully-hadronic decay 
channel arises from the hypothesis, used throughout the likelihood construction, that 
the fragmentation of the partons from both W bosons happens independently. However, 
Bose- Einstein Correlations (section 6.11) and colour reconnection (section 6.12) effects 
may result in cross-talk between the two W systems. A jet reconstruction technique is 
presented here which has been designed to have reduced sensitivity to colour reconnection 
effects. 

Conventionally, as used for the jets in the semi-leptonic analysis, the particles in the 
event are clustered into jets using a jet clustering algorithm and the energy, magnitude 
of the momentum and direction of the jet are reconstructed from the clustered particles. 
The jet momentum and energy are then used as the input to the kinematic fit. This 
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Figure 4: The left hand plots show the distribution of the Dp^r variable for four jet 
(top) and five jet (bottom) events from the full Lep2 data sample and the correspond- 
ing simulation samples. The right hand plots show the distribution of the four-fermion 
event purity with this variable at a centre-of-mass energy of 199.5 GeV extracted from 
simulation events. The fitted parameterisation of this distribution is given by the line. 
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Figure 5: Illustration of the iterative cone jet reconstruction algorithm used for the 
fuUy-hadronic W mass analysis as discussed in the text. 

technique is referred to in this paper as the standard reconstruction method and provides 
the optimal statistical sensitivity. 

In the alternative reconstruction algorithm discussed here the effect of particles in 
the inter-jet regions on the reconstructed jet direction is reduced. This is achieved by 
using a cone algorithm. The initial jet direction p^^* is defined by the standard clustering 
algorithms (DURHAM [27], CAMBRIDGE [28] or DICLUS [29]) and a cone of opening angle 
Rcone defined around this as in figure 5. The jet direction is recalculated (direction (1) 
on the figure) using those particles which lie inside the cone. This process is iterated 
by constructing a cone (of the same opening angle Rcone) around this new jet direction 
and the jet direction is recalculated again. The iteration is continued until a stable 
jet direction p^^^^ is found. Only the jet direction is changed in this procedure, the 
magnitude of the momentum and the jet energy are rescaled to compensate for the lost 
energy of particles outside the stable cone. The value of the cone opening angle Rcone is 
set to 0.5 rad, a value optimised for the measurement of the colour reconnection effect as 
reported in [30]. 

This cone jet reconstruction technique reduces the sensitivity to the colour reconnec- 
tion effect (see section 6.12) at the expense of some statistical sensitivity. The expected 
statistical uncertainty increases by approximately 14%. This technique has been applied 
only to the W mass and not to the W width analysis. 

This technique of jet reconstruction should not be confused with the alternative jet 
clustering algorithms (DURHAM, CAMBRIDGE or DICLUS) used in the analysis (see below). 
The alternative jet clustering algorithms are used as the starting point for the cone jet 
reconstruction in order to improve the statistical sensitivity of the analysis rather than 
to reduce the sensitivity to colour reconnection effects. 
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5.3.3 Likelihood Function 
Event Ideograms 

Each of the selected events is analysed through the use of a likelihood ratio function, 
which we will label here as the event ideogram. The final ideogram for each event consists 
of the weighted sum of the ideograms produced using a range of event reconstruction hy- 
potheses hi. These reconstruction hypotheses, including for example the possible different 
associations of the jets to their parent W bosons, are discussed below. The details of how 
these hypotheses are combined is then described below under the heading of 'Ideogram 
Sum'. 

The ideogram reflects the relative compatibihty of the kinematics of the event with the 
premise that two heavy objects, with masses nix and niy, were produced. The ideogram 

is based on the least-square, xlc^ °f energy and momentum constrained flt of the 
observed set of jet kinematics, {pj}, of the reconstructed final state. 

Thus, for each pair of test masses fh = {mx,my), we can obtain the xiciiPj}]''^^ ^i)- 
As the calculation of the over the full mass m plane is computationally intensive we 
apply the following approximation in the analysis. The is only calculated once per 
hypothesis hi at the minimum of the x|^(m) in the full m-space. The probability in all 
other points fn = {mx,my) is calculated using a Gaussian approximation for the x^l'^) 
given by: 

xUmx, my) ~ xlc + (m - - ra^'), 

with 

_ / nT,x \ 

The masses m^'*, m^**, and the covariance matrix V are taken from the 4C kine- 
matic flt. When the xlc larger than the number of degrees of freedom (NDF=4), the 
xi{mx,my) is rescaled with a factor NDF/xl^^ in order to compensate for non-Gaussian 
resolution effects. 

This procedure decreases the computing time taken by an order of magnitude com- 
pared with the full six constraint fit [3], while resulting in only a minimal reduction in 
the W mass precision obtained (2 ± 1%). 

We denote the ideogram of the event under hypothesis hi as P({pj}|m, hi). Assuming 
a Gaussian form, this is calculated from the x^ as follows: 

P{{pj}\m,hi) dm = e-Ki> (^-^ xlc{{Pj}\m,hi)^ dm. 

Example ideograms are shown in figure 6. These ideograms show the weighted sum of 
the reconstruction hypothesis ideogram terms for an individual event. The reconstruction 
hypotheses, which we will discuss in the following sections, include a range of options for 

the jet clustering algorithms that assign particles to jets, the possible associations of jets 
to W bosons, and a treatment for events that may have significant initial state radiation. 
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Figure 6: Examples of a reconstruction hypothesis weighted sum of two-dimensional prob- 
ability ideograms (see text) for a four-jet (left) and five-jet (right) hadronic event. The 
ideograms include terms from each potential jet-pairing, three jet clustering algorithms 
and possible ISR emission. The 1,2,3 and 4-sigma contours are shown. 



Jet Pairings 

As discussed in section 5.3.1, the reconstructed particles in the event were clustered 
into four or five jets. These jets can then be associated to their parent W bosons. For 
events clustered into four (five) jets there are three (ten) combinatorial possibilities for 
the jet pairing into W bosons. The relative probability of each of these jet pairings to be 
the correct jet association is estimated. 

This jet to W boson association weight, Wk-, is estimated as a function of the recon- 
structed polar angle of the W boson and the estimated charge difference between the two 
reconstructed W bosons in the event. For the five jet events the transverse momentum 
of the gluon jet is also used. 

The production angle 9w of the W"*" (W~) boson is correlated with the flight direction 
of the incoming (e~) beam. For each jet pairing the W boson polar angle was calculated 
and its probability PeiOw) assessed from a centre-of-mass dependent parameterisation of 
correctly paired simulation events. 

The jet charge Q*g^ for jet i in the clustered event can be measured as: 



Q]et — 
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where njet are all charged particles in jet i, while qn and pn are their charge and momen- 
tum. For each association k of the jets to their parent W bosons the charge difference 
'^Qk = Q'k^ " is obtained. Again, the probability of this being the correct jet 
assignment is assessed using a Monte Carlo simulation-derived parameterisation. The 
relative weight for each jet pairing k can be expressed as: 

= Pw^A Qk) ■ Pe(4j + (1 - Pw+{A Qk)) ■ Pei-^ - 

In five jet events, a two jet and a three jet system are considered. The three jet system 
is considered as comprising a qq pair and a gluon jet. The probability of emission of a 
gluon from a qq pair is approximately inversely proportional to the transverse momentum 
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of the gluon with respect to the original quarks. Hence, the most probable gluon jet in 
the three jet system is the jet with the smallest transverse momentum (kx) with respect 
to the two other jets in the candidate W boson rest frame. Each of the ten possible jet 
associations, in this five jet event, is then given a relative weight from its most probable 
gluon jet of wl''^°^ = l/kx- 

The combined relative jet pairing weight of each combination is given by multiplying 
the jet pairing weights and, for five jet events, also multiplying by the wl''^""' weight. 
The relative weights are then normalised so that the sum of the weights for all the jet- 
paring combinations of the event is 1, giving combination weights Wk- The use of all the 
jet pairings, rather than simply picking the best one, improves the statistical precision of 
this analysis by 4%. 

Jet Clustering Algorithms 

Several standard jet clustering algorithms are used in this analysis. Whilst the overall 
performances of the algorithms are similar, the reconstruction of an individual event can 
differ significantly. In this analysis, the event ideograms were reconstructed with three 
clustering algorithms DURHAM, CAMBRIDGE and DICLUS. The ideograms resulting from each 
clustering algorithm arc summed with fixed optimised relative weights, Wc- determined 
from simulation events. The sum of the three jet clustering weights for one event is 
normalised to 1. 

The use of a range of jet clustering algorithms, rather than taking only one, improves 
the statistical precision of this analysis by 5%. 

Initial State Radiation Hypotheses 

A kinematic fit (see section 5.1) is performed with modified constraints and an extra 
free parameter pf* to account for the possible emission of an ISR photon of momentum 
Pz inside the beam pipe. The modified constraints are: 

'^objects 

1=1 

The probability that the missing momentum in the z direction is indeed due to an 
unseen ISR photon was extracted from the simulation as a function of where 
(7p^ is the estimated error on the fitted z momentum component; only events with this 
ratio greater than 1.5 are treated with the mechanism described below. 

Additional ideograms are then calculated for these events, with a relative weight factor 
derived from the ISR hypothesis probability. The ideogram obtained without the ISR 
hypothesis is given a relative weight 1, while the other ideograms obtained from this 
procedure are given relative weight factors according to the distribution shown in figure 7. 
The weights are then normalised such that the sum of the ISR and no ISR hypotheses 
for an event sum to 1, giving ISR weights Wisr- 

This treatment is applied to 15% of the events and results in an improvement of the 
expected W mass error for these events of 15%. 

Ideogram Sum 

An ideogram is produced for each event under each of the possible reconstruction 
hypotheses. For four jet events there arc three jet association hypotheses to be performed 
with three clustering algorithms and maximally two ISR hypotheses, giving a total of 
eighteen ideograms. For five jet events there are sixty possible ideograms. The final 
ideogram for each event is produced as a weighted sum of these: 
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Figure 7: Parameterised weight given to the ISR solution of the kinematic fit, relative 
to the unity weight of the no ISR solution, as a function of the |pf*|/crp^ value of the 
event for different centre-of-mass energies. The period with a damaged TPC sector (S6) 
is indicated with a dashed line. 



3 or 10 2 3 

P{{pj}\m, {hi}) = J2 J2 ■ Wisr ■ Wc- P{{pj}\m,hk,isr,c), 

k=l isr=l c=l 

where the sum over k takes into account the three or ten possible jet pairings in the event, 
the sum over isr the two different initial state radiation hypotheses used in the kinematic 
fit and the sum over c the three jet clustering algorithms. The sum of all weights for 
each event is fixed to unity, so that while possible reconstruction hypotheses within an 
individual event have different weights the overall weight for each event is the same. 

Likelihood 

To obtain information about Mw and Fw a theoretical probability distribution func- 
tion, P(m|Mw,Fw), is required predicting the population density in the m-plane of the 
event ideogram. The ideogram in m-space can then be transformed into a likelihood, 
£e(Mw,Fw), in the (Mw,Fw)-space by convoluting it with this expected distribution 
P(m|Mw,Fw): 

£e(Mw,Fw)=/ / P(fe}|m,{/i,})-P(m|Mw,Fw) rfm, (6) 
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where the two-dimensional integral is over the relevant kinematic region in the m-space. 
This region is taken to be mmin = 60 GeV/ and mmax = 110 GeV/ c^, and the combined 
ideogram is normalized to unity in the same region: 

/ / P{{pj}\m,{hi}) dm = 1. 
Theoretical Distribution Function 

The theoretical probability distribution function, P(m|M-w,rw), predicts the popula- 
tion density in the m-plane of the event ideogram for a given Mw and Fw To provide 
an accurate description of the data the form assumed for P(m|Mw, Tw) must take into 
account not only the expected distribution for the W+W" — > qq'qq' signal events but 
also that of the background events in the selected sample. The two principal components 
of the background, Z — > qq{j) and ZZ qq'qq', are considered. 

The background process Z — qq{'j) does not have a doubly resonant structure and 
a uniform population of these events is expected in the m-space independent of the 
values of the parameters (Mw,rw)- Therefore, the probability density function from this 
background source is assumed to be a constant denoted B. The probability (P^-^) that 
a given event is a qq'qq' event was calculated from the event topology as described in 
section 5.3.1. 

The W+W^ qq'qq' and ZZ — * qq'qq' events both have a doubly resonant Breit- 
Wigner structure in the m-plane, modulated by a phase-space correction factor PS{m\^/s) 
due to the nearby kinematic limit mw+ + mw- < V^- The probability density function 
component used to model four-fermion events is given by: 

5(m|Mw,rw) = P5(m|V^)- 

~ww ~zz 

• -/vv+^zz • PWVw(m|Mw, Tw) + ^wwV^zz • PW^zz(m|Mz, Tz)] , 

where and reflect the accepted cross-sections, calculated from simulation, of 

respectively the W+W^ and the ZZ final states. These cross-sections are centre-of-mass 
energy dependent but arc independent of the reconstructed event topology. 

The two-dimensional Breit-Wigner distribution is approximated as the product of two 
one-dimensional Breit-Wigners: 

5WW(m|Mw,rw) ^ BWw{mw+\Mw,r^) ■ BWw{mw-\My^,r^), 

with BWy/ given by the expression in equation 5 of section 5.2.3. An expression of the 
same form is assumed for the ZZ component. 

A dependence on the centre-of-mass energy is also introduced into S'(m|Mw,rw) 
through the phase space correction factor PS{m\^/s)■. 

PS{m\^) = ^y'(s-m^+ - m^-Y - Am^+m^-. 

The combined density function is then constructed from the signal and background 
terms: 

P(m|Mw, Tw, v^) = P^^ • -S(m|Mw, Tw, ^i) + (1 - P^O • B. 

Utilising this probability density function, and the event ideogram, equation 6 may 
be used to calculate the event likelihood function. The extraction of the parameters of 
interest, Mw and Fw, from the event likelihood functions are discussed below. 
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5.4 Mass and Width Extraction 

The mass and width of the W boson are extracted from maximum hkehhood fits to 
data samples. This section describes this procedure, the cahbration apphed and the 
cross-checks of this method that have been performed. 

The distribution of the reconstructed invariant masses of the selected events after 
applying a kinematic fit, imposing four-momentum conservation and the equality of the 
two di-jet masses, are shown in figure 8. This figure is provided for illustrative purposes 
only, the mass and width fitting procedure is described below. 

The combined likehhood of the data can be obtained from the product of the event 
likelihoods described above. In practice this is achieved by performing the sum of the 
logarithms of the individual event likelihoods. The fitted data samples are divided by 
data taking year and applied event selection. For the mass fit the data from the fully- 
hadronic event selection and the electron, muon and tau semi-leptonic selections are all 
fitted separately. In the determination of the W width, where the relative precision is 
much worse, the data arc divided only into fully-hadronic and semi-leptonic selection 
samples. The procedure for combining the results from each of these fits is discussed in 
section 7. 

The W mass and width are extracted from maximum likehhood fits. The W mass fit 
is performed assuming the Standard Model value for the W width (2.11 GeV/c^). The 

W width was obtained assuming a mass of 80.4 GeV/c^. The correlation between Mw 
and Fw was found to have a negligible impact on the extracted mass and width value: 
the current uncertainty of 44 MeV/ on Fw [32] gives rise to a 0.6 MeV/ uncertainty 
in the extracted Mw- 

The terms used in the likelihood and described above are functions which approximate 
a description of the underlying physics and detector response. Hence, this approach 
necessitates a calibration of the analysis procedure. The calibration is performed using 
signal and background simulation events for which the true mass and width values are 
known. Rather than regenerating the events at a range of mass and width values, the 
calibration of the analysis uses reweighted events. The reweighting was performed using 
the extracted matrix clement of the WPHACT and YFSWW generators. The reweighting 
procedure is cross-checked using independent simulation events generated at three W 
mass and width values. In the fully-hadronic channel where both the standard method 
and the cone-jet reconstruction technique are applied to the W Mass measurement, both 
analyses are calibrated separately: the illustrative values reported in this section are for 
the standard analysis. 

A high statistics simulation sample is used to calibrate the analysis, comprised of an 
appropriate mixture of signal and background events. The result of the likelihood fit as 
a function of the simulated W mass is shown in figure 9 for the yuz7^qq' channel analysis 
at ^/s = 189 GeV. The analysis has a linear behaviour in the mass window of interest, 
and the calibration curves are defined by two parameters : 

• the slope of the generated mass against fitted mass line; 

• the offset defined at a fixed reference point. This point is chosen to be the value 
used in our simulation; 80.4 GeV/ for the mass and 2.11 GeV/ for the width. 

The slopes at different energies are found to be compatible, and their mean values 
are respectively 0.984 ± 0.013, 0.993 ± 0.006 and 0.963 ± 0.013 in the eV^qc^, /xl^^qq' and 
TVrq^ analyses. In the qq'qq' analysis the slope was compatible with unity to within 2% 
at all centre-of-mass energies and no slope calibration was applied. 
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Figure 8: The distribution of the reconstructed W masses from a kinematic fit with 
five constraints imposed in the (a) el/eqq', (b) /iF^qq', (c) TVT-Cf( and (d) and (e) qq'qq' 
analysis channels at all energies, (d) shows the data sample taken at all energies until 
September 2000, the data taken after that with a damaged TPC sector is shown in (e). 
In (d) and (e) only the jet pairing with the highest probability is included in the figures. 
The simulation samples have been normalised to the same integrated luminosity as the 
data. 
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Figure 9: W mass calibration curve in the /xz^^^qq' channel at = 189 GeV. The dashed 
line indicates the result that would be obtained without any analysis bias. 

The highly linear behavior, with a value of the slope close to unity is an a posteriori 
justification of the fitting functions used in the likelihood fit and described in section 5.3.3. 
The remaining effects not taken into account by these fitting functions give rise to the 
offset. As an example, the calibration offsets at ^/s = 189 GcV are respectively — 0.108 ± 
0.012, -0.215 ± 0.010, -0.252 ± 0.015 and -0.222 ± 0.006 GeV/c^ in the eU^q^, iiV^q^, 
Ti^T-qq' and qq'qq' analyses for the mass. The offsets vary slightly with the centre-of-mass 
energy. 

The same procedure is also applied for the W width analyses. In the iu^qq' channel 
a slope of 0.894 ± 0.008 is obtained independent of the centre-of-mass energy and the 
offset at ^/s — 189 GeV was -1-0.065 ±0.015 GeV/ . However, in the qq'qq' analysis the 
slope is found to be dependent on the centre-of-mass energy, the slopes at ^/s — 189 GeV 
and 205 GcV are approximately 1.1 and 1.2 respectively and furthermore the relation 
between the reconstructed and generated Fw is not perfectly linear. Hence the offset 
is parameterised as a function of the generated W width and the centre-of-mass energy. 
The cahbration offset at ^/s ^ 189 GeV is 183 ± 13 MeV/c^ at the reference width. 

The analyses are corrected with these calibration results, and the statistical error on 
the offset is included in the systematic error (sec below). 

After applying the calibration procedure, the consistency of the analyses is checked. 
Sets of simulation events, with a sample size the same as the data, containing the expected 
mixture of signal and background events were used to test the analyses. Figure 10 shows 
error and pull plots from analysing 20000 or more such samples, where the pull is defined 
as 

_ (Mwfit ~ Mwgen) 

here the subscript 'fit' and 'gen' distinguish the result from the calibrated analysis fit and 
the generated parameter in the simulation respectively. The cxfit is the error estimated 
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by the analysis. This error has been scaled in the analysis to obtain a Gaussian width 
of one for the pull distributions, as shown in the plots. These plots were produced at 
all centre-of-mass energies for both parameters. The error distributions in figure 10 also 
demonstrate that this quantity is in good agreement with the value obtained from the 
data. 
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Figure 10: The errors (left) and pulls (right) of the W mass fits for each semi-leptonic 
analysis channel and the fuUy-hadronic channel. These plots were obtained using simu- 
lated event samples with the same statistics as the data sample collected at 200 GeV. The 
errors obtained on the fits to the data samples were 365 MeV/ (? for the eX/gqq' analysis, 
282 MeV/c^ for /xI7^qq', 438 MeV/c^ for rl/^qq' and 149 MeV/c^ for the standard qq'qq' 
analysis. 
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6 Systematic Uncertainties 

The sources of systematic error that have been considered for the W mass and width 
determinations are described in the subsections below. The results of these studies at 
example centre-of-mass energies are summarised in tables 14, 15 and 16. In the fully- 
hadronic channel the standard method and the cone jet reconstruction technique have 
been applied as described in section 5.3.2. The systematic uncertainties are in agreement 
between these two techniques except for the error sources from final state interactions 
(FSI), where separate values for the two techniques are given. 



6.1 Calibration 

The analysis calibration procedure is described above in section 5.4. The accuracy 
with which the offset of the analyses can be determined is limited by the size of the 
generated simulation samples. Sufficient events were generated to limit this error to 5% 
or less of the statistical error on the mass or width determination in any given channel. 



6.2 Detector Effects - Muons 

Contributions to the systematic error on the W mass and width due to the reconstruc- 
tion of muons are considered in this section. These were evaluated using the Z — > /i"^//" 
events collected at the Z peak during the Lep2 period. The systematic uncertainties 
determined by these studies for the W mass analysis are presented in table 3. 

Inverse Momentum Scale 

The primary sources of systematic error on the muon momentum scale are the detector 
alignment or possible reconstruction distortions (particularly in the TPC). As a result of 
these effects, we may also anticipate an opposite bias on the measured track curvature 
for positive and negative muons. 

Corrections to the inverse momentum scale, l/p, are calculated from the selected /x^/x" 
samples. The mean inverse momentum, < 1/p >, is calculated separately for positive 
and negative muons in different bins of the polar angle, and a correction for the positive 
muons is defined as 

^(< — > - < 4: >)' (7) 

2 p— p+ 

with the opposite sign correction applied to negative muons. These corrections are typi- 
cally of the order 1 to 2 x 10"*^ GeV~^c, except in the polar angle regions at the junction 
between the barrel and endcaps where the correction can reach 10~^ GeV~^c in the worst 
case. In the simulation this correction is, as expected, compatible with zero. After apply- 
ing the corrections < 1/p >data and < 1/p > simulation are found to be in agreement within 
0.2%, and this value is used to calculate the systematic on the muon inverse momentum 
scale. The systematic uncertainty on the positive and negative muon inverse momentum 
scale difference is estimated by varying the correction by ±50% of its value. 

Inverse Momentum Resolution 

The momentum resolution (typically 0.001 GeV~^c in was found to be commonly 
around 10% better in simulation events than in the data. This discrepancy, determined 
for all years of Lep2 and polar angle regions, is corrected by smearing the simulation with 
a Gaussian. An additional smearing of ±0.0003 GeV~^c in 1/p is used to estimate the 
systematic error resulting from this correction. This systematic does not affect the Mw 
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determination but is a small component of the Fw measurement uncertainty for events 
containing muons. 

6.3 Detector Effects - Electrons 

Contributions to the systematic error on the W mass and width due to the reconstruc- 
tion of electrons are considered in this section. These were evaluated using the Bhabha 
and Compton events collected at the Z peak and high energies during the Lep2 period. 
The systematic uncertainties determined by these studies for the W mass analysis are 
presented in table 3. 

Energy Scale 




DELPHI 



.f f 




1 

E/Ebeam 



Figure 11: The ratio E/Ebeam for electrons in the endcaps from Bhabha events recorded at 
the Z peak in 1998. The shaded histogram is the simulation and the points are the data. 
Plot (a) shows the raw distribution, while plot (b) gives this after the bremsstrahlung 
correction discussed in the text. The resolution correction (see text) has also been applied. 



The reconstructed energy of electrons was compared between data collected at the Z 
peak and fully simulated samples of Bhabha events. In the barrel region of the detec- 
tor the data and simulation are in good agreement. However, in the forward directions 
a slight difference is observed between the data and simulation (see figure 11) and at- 
tributed to an under-estimation of the quantity of material in the simulation before the 
electromagnetic calorimeter in the Delphi endcaps. A correction is applied to the sim- 
ulation by introducing the effect of extra bremsstrahlung emission corresponding to an 
additional 3% of a radiation length. Following [33], the probability w that an electron of 
initial energy Eq has an observed energy between E and E+dE after traversing a thickness 
of t radiation lengths is 

u,(^„.^..).^^f '""X)'7- . (8) 

Eq T{t/ln2) 

For each event, the corrected energy E is chosen randomly according to the distribution 
w. The optimal value of the parameter t was adjusted from the data and simulation 
comparison. 
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After the endcap correction was applied, good agreement between data and simulation 
was obtained throughout the detector. The residual systematic error on this absolute 
energy scale is estimated to be ±0.3% of the measured energy and is estimated from the 
selection cut stability and statistical precision of the data and simulation comparison. 

Energy Resolution 

The resolution on the reconstructed electron energies was also compared between the 
data and simulation Bhabha samples. The agreement is improved by applying a Gaussian 
smearing to the simulation with a width varying between 1 and 2% of the measured 
electron energy in the barrel, and 2 to 4% in the endcaps, depending on the year of data 
taking. The systematic error on this smearing Gaussian width is estimated to be ±1% 
of the measured energy. This systematic does not affect the Mw determination but is a 
small component of the Fw measurement uncertainty for events containing electrons. 

Energy Linearity 
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Figure 12: The double ratio of reconstructed and true average energy values in data and 
simulation, ( Eyec/Etrue )data/{ Erec/ Etrue )mci fo^ data taken in 2000. The shaded area 
represents the quoted systematics due to a possible dependence of the energy calibration 
with the electron energy. The left hand plot is for electrons observed in the barrel 
electromagnetic calorimeter and the right hand plot for electrons in the endcap. Note 
that, by construction, the Bhabha point at 45 GeV is at one. 

The reconstructed electron energy was also studied as a function of the true energy. 
The Z peak and high energy running provided high statistic Bhabha samples with which 
to study electrons of 45 GeV and above 100 GeV energy. For these samples the "true" 
electron energy is taken from the beam energy. The reconstructed electron energy was also 
checked using low energy electrons from Compton events at the Z peak, and high energy 
electrons from radiative Bhabha scattering at high centre-of-mass energy. In these cases 
the true energy of the lepton is deduced from 3-body kinematics using only the angular 
information and assuming that the unseen particle was along the beam axis. Figure 12 
shows the compatibility of the reconstructed electron energy in data and simulation, only 
statistical errors are shown. One of the three points measured for radiative Bhabhas in 
the Barrel shows a discrepancy but this effect is not confirmed by the better measured 
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high energy (non-radiative) Bhabha point, whereas physical cahbration problems such as 
threshold effects or leakage in the calorimeter would be expected to increase in size with 
energy. Hence, no additional corrections are applied. A systematic error is estimated 
assuming a deviation of the energy calibration slope Edata/Egimuiation versus Esimuiation of 
1% over the range 25 to 70 GeV. These values approximately correspond to the relevant 
energy range for the observed electrons in the analysis. 

6.4 Detector Effects - Taus 

The rl/^qq' channel differs from the other W"'"W~ semi-leptonic decay channels as 
these events contain two (or three for leptonic tau decays) neutrinos in the final state. 
Thus, the mass of the event can be determined only from the decay products of the other 
W. As a result the lepton systematics described in the preceding sections are not relevant 
to the rl/T-qq' channel. The only relevant systematic involving the tau decay products 
arises from uncertainties in the assignment of the reconstructed tracks between the tau 
product and the hadronically decaying W. This effect is small compared with the overall 
uncertainty on the jet energy and direction, the systematic on which is considered in the 
sections below. 



Mw Lepton Correction Systematic Errors (MeV/c^) 


Sources of Systematic Error 


cI7eqq' 189 GeV 


el/^qq' 205 GeV 


Electron Energy Scale 
Electron Energy Resolution 
Electron Energy Linearity 


18 

16 


22 

11 




//I7^qq' 189 GeV 


/x77^qq' 205 GeV 


Muon 1/p Scale 

fx^ fx^ 1/p Difference 

Muon 1/p Resolution 


16 
1 


21 
4 



Table 3: Contributions to the systematic error on the W mass measurement at 189 and 
205 GeV related to the lepton reconstruction. The uncertainties on each of these numbers 
is typically 3 MeV/c^. 



6.5 Jet Description 

Jets are composite objects, and the detector and analysis response to them can be 
dependent on their internal structure. Therefore it is not straightforward to separate in a 
clean way uncertainties arising from the modelling of the detector in the simulation from 
those due to the theoretical description of the jet structure. 

Moreover this description is not based on exact calculations, whose uncertainty can 
be in principle reasonably well estimated, but on phenomenological models tuned to best 
reproduce the data at the Z peak: the Lund model as implemented in PYTHIA is the 
standard choice for this analysis. In this situation the comparison of different models 
may be a useful tool to understand which parts of the fragmentation description the 
measurement is sensitive to, but only a direct comparison of the chosen model with well 
understood data samples, in particular Z hadronic decays, can give the ultimate estimate 
of the uncertainty from the observed data-simulation disagreements. 
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The jet studies performed are described in the text below and the corresponding jet 
correction systematic errors are provided in table 4. The most relevant jet characteristics 
were calibrated on real data control samples, and uncertainties on these calibrations are 
propagated through the analysis. 

Energy Scale 

The absolute jet energy scale was studied in on-peak Z — > gg decays, by comparing 

the reconstructed energies, Erec, in data and simulation in selected two jets events. The 
6— tagging technique is used to remove b quark jets which are essentially not present in 
W+W^ decays. The true jet energy in these events is assumed to be the beam energy 
Eheam-i Under the assumption that the bias introduced by QED ISR is described with 
negligible error in the simulation (the KK2f generator was used for these events). The 
double ratio of average values ( Erec/Ei,eam )data/{ Ej.ec/Ebeam )mc was evaluated as 
a function of the jet polar angle and applied as a scale factor correction to the four- 
momentum components of the jet in simulated events. The correction value depends on 
the year as well as the angular region, with the deviation from unity ranging typically 
from a few per mille up to 3-4% in the most forward region. 

The systematic uncertainty on this correction is determined by the limited on-peak Z 
statistics, and it is estimated to be ±0.3%. 

Energy Resolution 

The same event sample used to study the jet absolute energy scale was also used to 
calibrate the jet energy resolution in the simulation. A Gaussian smearing was determined 
from the data and is applied to the simulated jet energy with a magnitude dependent on 
the ratio of the reconstructed and true jet energies. This procedure takes into account 
the asymmetric shape of the jet energy observable. When applying the correction to the 
simulated W"'"W~ events an estimate of the true jet energy is required. When the event 
is reconstructed with two jets from each hadronically decaying W, the generated quark 
energies are used. However, when gluon radiation has given rise to an additional jet the 
true jet energy estimate is determined by applying the same clustering algorithm as used 
in the analysis to the simulated partons prior to the detector simulation. In both cases 
the association of the true and reconstructed jets is performed according to geometric 
criteria. 

The average resolution correction ranges from 4.5% of the jet energy in the barrel 
to 6.6% in the cndcaps. The correction is also dependent on the year. The systematic 
uncertainty on the correction is estimated to be ±2% of the jet energy. 

Energy Linearity 

The dependence of the energy calibration as a function of the jet energy was checked 
using low energy jets from qq -\- gluon events at the Z peak and high energy jets from 
e^e" — > qq decays at high energy. 

In the first case, the true jet energy is determined Tising three-body massless kinemat- 
ics. The jet energy range used in this study is restricted to the region where the data 
and simulation true energy distributions do not show sizeable discrepancies. This energy 
selection avoids introducing an unnecessary sensitivity in this analysis to the modelling 
of hard gluon radiation in the simulation. 

In the second high-energy jet case the effective hadronic mass is required to be 
such that \J s' / s > 0.95. The true jet energy is then again determined using three-body 
massless kinematics but now the third object is an hypothetical ISR photon emitted along 
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the beam pipe. The difference between the estimated jet energy and the nominal beam 
energy is constrained to be smaller than 10 GeV. 

A jet energy linearity slope in Edata/Egimuiation versus Egimuiation is then determined. 
The study was performed separately in the barrel and cndcap regions of the detector and 
for each data taking year. The results from the different data taking years are compatible 
within statistical errors. The study showed agreement in the slope at typically the 0.5% 
level over the range 25 to 75 GeV, and this deviation value is used to determine the 
systematic uncertainty. 

Angulcir BicLS 

As reported in [9] , the reprocessing of data and simulation used for this analysis has a 
noticeable excess of tracks at low polar angles (forward tracks) in data as compared to the 
simulation. The most likely cause of this effect is an underestimation in the simulation 
of the track reconstruction efficiency for low-momentum particles at low polar angle. 

This effect introduces a small bias in the distribution of the jets' reconstructed polar 
angle in the simulation compared with data. In order to evaluate the effect of such a 
bias, a systematic shift of the jets' polar angle is applied to the simulation. The shift as 
a function of the polar angle itself has been determined using on-resonancc Z hadronic 
decays, and is found to have the form 0.008 cos where < 6'j < 7r/2 is the polar 
angle of the jet. The corresponding W mass and width shifts have been evaluated and 
symmetric systematic errors of these values applied. The W mass uncertainty is reported 
in table 4. 

Anguletr Resolution 

A study of the acoUinearity of jets in on-peak Z ^ qq events was performed and 
appropriate smearings to the simulation of the jet angular direction, dependent on the 
polar angle of the jet, were estimated. The smearings on the polar angle are typically 
5 mrads. A systematic error is estimated by applying an extra 5 mrad angular smearing. 

Jet Mass 

The jet mass is known not to be exactly described in the simulation; both inaccuracies 
in the fragmentation description (related to the jet breadth due to soft and hard gluon 
radiation) and imperfections in the modelling of the detector response (reconstruction 
efficiencies and noise) are responsible for these discrepancies. However, only those data- 
simulation differences in the jet mass which are not compensated by differences in the 
inter-jet angle are relevant for the systematic uncertainty, since these cause systematic 
biases in the reconstructed W mass. 

For this reason the fragmentation-induced differences are only marginally relevant for 
the mass measurement. Furthermore, the calibration procedure adopted, in particular for 
the energy and angular smearing, corrects for most of the effects given by the differences 
in jet breadth. The jet breadth is relevant as broader jets are worse reconstructed: they 
are detected with larger uncertainties on the jet direction; are likely to lose more energy 
due to the imperfect hermeticity of the detector; and cause more confusion in the jet 
clustering. 

The jet correction procedure described above, as well as the constrained kinematic 
fit, modifies all the four-momentum components of the jet but leaves unchanged the jet 
boost, i.e. the E/m ratio. It is therefore useful to study this observable, instead of the 
simple jet mass. 

Detector noise is a source of data-simulation discrepancy which clearly biases the 
reconstructed boson mass, since it changes the mass and boost of the jets while leav- 
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ing, on average, the inter-jet angle unchanged. Significant data-simulation differences 
in low energy neutral clusters, both in the electromagnetic and hadronic calorimeters, 
are attributed primarily to an imperfect noise description, while the discrepancies in the 
charged particles of jets are considered to be almost entirely due to the modelling of the 

fragmentation. 

The average effect of removing low energy neutrals below 2 GeV on the jet m/E 
was evaluated as a function of the polar angle and of the m/E ol the jet itself, since the 
impact of the noise depends on the breadth of the jet. The expected effect on the neutrals 
from fragmentation was subtracted. The fragmentation effect was obtained from charged 
particles, suitably scaled for the relative neutral and charged particle multiplicity. 

This m/E effect was then propagated in the full analysis chain to extract the relative 
systematic uncertainty on the full mass and width measurements. 



Mw Jet Correction Systematic Errors (MeV/c^) 



Sources of Systematic Error 


189 GeV 




ez/eqq' 


/iz/^qq' 


rz/^qq' 


qq'qq' 


Energy Scale 


8 


6 


11 


8 


Energy Resolution 


3 


3 


5 


9 


Energy Linearity 


12 


9 


12 


16 


Angular Bias 


3 


5 


5 


2 


Angular Resolution 








8 


Jet Mass 


9 


8 


8 


10 






205 GeV 






ei/eqq' 




Ti/^qq' 


qq'qq' 


Energy Scale 


11 


9 


16 


8 


Energy Resolution 


8 


5 


8 


10 


Energy Linearity 


15 


11 


20 


8 


Angular Bias 


9 


8 


7 


19 


Angular Resolution 








1 


Jet Mass 


13 


12 


17 


13 



Table 4: Contributions to the systematic error on the W mass measurement at 189 and 
205 GeV related to jet reconstruction. The uncertainties on each of these numbers is 
typically 6 MeV/c^. 



Fragmentation Model 

The effect of using different hadronisation models on the analysis was studied by 
replacing the standard choice, PYTHIA, with both the ARIADNE and HERWIG models, each 
tuned by DELPHI to best match experimental data. The mass and width shifts were 
evaluated at 189 GeV and 207 GeV centre-of-mass energies and are reported in tables 5 
and 6. Detailed studies performed at the Z peak showed that for several observables all 
the models showed disagreements with the data and that these disagreements were all 
in the same direction: the jet mass variable, discussed in the previous paragraph, is a 
clear example. Hence the results of the hadronisation model comparison were used only 
to investigate the sensitivity of the analysis to specific features of the models, and not 
used directly as an evaluation of the systematic uncertainty due to the choice of model. 
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eI7eqq' 


AMw 
/iI7^qq' 


MeV/c^ 
rI7^qq' 


qq'qq' 


HERWIG - PYTHIA 
ARIADNE - PYTHIA 


-7± 10 
-11 ±9 


-16 ±9 
-12 ±9 


-17± 13 
-10± 12 


-9 ±5 
-15 ±5 



Table 5: Effect of different fragmentation models on the W mass determination. 





AFw M 
£z7£qq' 


[eV/c" 
qq'qq' 


HERWIG - PYTHIA 
ARIADNE - PYTHIA 


+46 ± 13 
-9± 15 


-2± 11 
+1± 11 



Table 6: Effect of different fragmentation models on the W width determination. 



The biggest difference was found to be between PYTHIA and HERWIG, and was shown to 
be largely due to the different production rates of heavy particles, mainly kaons, protons 
and neutrons. At parton level these differences modify not only the jet masses but also 
change the jet-jet angles accordingly, leaving the bosons invariant masses unchanged. 
However, the reconstruction and analysis procedure breaks this compensation since in 
the fully-hadronic event reconstruction all charged particle tracks are assigned the pion 
mass, and all neutrals are assumed to be massless (photon- like) . In the semi-leptonic 
analysis, the nominal masses are used in the jet reconstruction for those particles with a 
positive identification, i.e. for charged kaons and protons identified by the RICH and for 
and Lambdas reconstructed as secondary vertexes from their decay products [10]. 

The HERWIG version used, although tuned to best reproduce the Z peak DELPHI data, 
is known to describe the particle production rates poorly. This is especially the case for 
baryons, therefore using HERWIG accentuates this particle mass assignment effect. Gener- 
ally the measured particle rates arc closer to those in PYTHIA and ARIADNE. Reweighting 
in the models the production rates of the most abundant heavy particles species, kaons 
and protons, reduces the disagreement among the different models, bringing it to the level 
of the statistical uncertainty of the fit. Tables 7 and 8 show the residual discrepancies 
obtained between the models after they have been reweighted to the PYTHIA values. The 
component of the fragmentation systematic error which is not due to the heavy particle 
multiplicity effect is obtained from these numbers. The largest vahic - either the central 
value or its uncertainty - from either model is taken as the systematic error estimate. 

The component of the fragmentation error that is due to the heavy particle rate was 
also evaluated for the W mass analysis; this small component of the error is neglected 
for the W width analysis. The W mass shift was evaluated between the DELPHI tune 
of PYTHIA and the same events reweighting to the measured particle rates ilcr of their 
uncertainty. The average of the modulus of the two shifts is reported in table 9 and is 
taken as the estimate of the fragmentation error due to the heavy particle multiplicity. 

The combined fragmentation error was evaluated for the W mass by adding the particle 
reweighting effects and the model variation uncertainty in quadrature. This fragmenta- 
tion error is listed separately from the other jet description uncertainties in the systematic 
uncertainty summary tables 14, 15 and 16. 
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AMw 
;Ul7^qq' 


MeV/c^ 


qq'qq' 


HERWIG Rew. - PYTHIA 
ARIADNE Rew. - PYTHIA 


-2 ± 10 
-10 ±9 


-8 ±9 
-10 ±9 


-5 ±13 
-10± 12 


-11±6 
-1±4 



Table 7: Effect of different fragmentation models on the W mass determination, after 
reweighting the heavy particle species rates in the Monte Carlo simulations to the mea- 
sured rates. 





AFw M 


eV/c" 
qq'qq' 


HERWIG Rew. - PYTHIA Rew. 
ARIADNE Rew. - PYTHIA Rew. 


+29 ± 13 
-11 ± 15 


+3 ±8 
-1±8 



Table 8: Effect of different fragmentation models on the W width determination, af- 
ter reweighting the heavy particle species rates in the Monte Carlo simulations to the 
measured rates. 

6.6 Mixed Lorentz Boosted Zs 

An alternative method of evaluating the jet description systematic is to use the tech- 
nique of mixed Lorentz boosted Zs (MLBZ). This method attempts to emulate W'^W" 
events using two on-peak Z events. The emulated W+W" events are constructed both 
from simulated events and the large statistics sample of Z peak data events. Standard 
W mass and W width analyses can then be performed on these event samples. Hence, 
the MLBZ method provides a direct comparison between data and the simulation model 
of choice. The difference between the measurements made from the data and simulation 
MLBZs can be interpreted as primarily providing a statistically sensitive cross-check of 
the fragmentation systematic assigned to the W mass and width measurements. This 
method would also identify some sources of detector modelling error. 







AMw MeV/c' 




Particle Type 


ez7cqq' 


/iz/^qq' 


rI7^qq' 


qq'qq' 




0.1 ±0.3 


0.9 ±0.3 


1.5 ±0.4 


0.2 ±0.5 


Proton 


2.0 ±0.4 


1.5 ±0.3 


3.2 ±0.5 


3.5 ±0.5 



Table 9: Effect on the W mass of reweighting the heavy particle species rates in the 
Monte Carlo simulations. The mass shifts were evaluated between the DELPHI tune 
of PYTHIA and versions reweighted to 1 sigma above and below the measured particle 
rates. The shift value reported is the average of the modulus of these two shifts. The 
measured charged multiplicity in a Z peak event for kaons is 2.242 ± 0.063 [32], whereas 
for protons the measured multiplicity is 1.048 ± 0.045 [32]. 
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A W"'"W~ — > ff 'f f event is emulated by selecting two Z events and rotating and Lorentz 
boosting them so that their superposition reflects a true W+W" event. The mixture of 
quark species will not be the same as in true W'^W" events, it will however be the 
same between the data and simulated Z samples that are used in the comparison. To 
emulate a qq'qq' event two hadronically decaying Z events were used. To emulate a iViqc^ 
event one Z decaying into hadrons and one Z decaying into charged leptons was used. 
One hemisphere of the Z— decay is removed to represent the W— >£|7^ decay. The 
emulation process is performed by manipulating the reconstructed tracks and calorimeter 
energy clusters. 

A realistic distribution of W"'"W~ events is obtained by using event templates. The 
four momenta of the four primary fermions in a WPHACT W+W" event are used as the 
event template. The Z events are chosen such that they have a thrust axis direction 
close to the polar angle of one of the W fermions. This ensures that the distribution 
in the detector of the tracks and energy clusters selected in the Z event follows that 
expected in W+W" events. Each of the template Ws is then boosted to its rest frame. 
The particles in a final state of a selected Z event are rotated to match the rest-frame 
direction of the fermions from the template W. The energy and momentum of the Z 
events arc then rcscalcd to match the kinematic properties of the W boson decay. The 
two Z events are then each boosted into the lab frame of the template W"*"W~ event 
and mixed together. The same W+W~ event templates are used for the construction of 
both the data and Monte Carlo simulation MLBZ events, thus increasing the correlation 
between both emulated samples. 

Tests were performed to confirm the reliability of the MLBZ method in assessing sys- 
tematic errors. MLBZs were produced using Zs with the PYTHIA, HERWIG and ARIADNE 
models and the observed mass shifts were compared and found to agree with the sta- 
tistically limited mass shifts observed in W''"W~ simulation events. A significant mass 
shift (300 MeV/c^) was introduced by using the cone rejection algorithm (discussed in 
section 5.3.2) for the W mass measurement in the qq'qq' channel. The real and simulated 
MLBZs and W+W" events agreed on the estimated size of the mass shift between the 
standard and cone estimators at the 15% level. 

The MLBZ method was used to create emulated W"'"W~ event samples. The Z events 
were selected from data recorded during the Lep2 calibration runs of the same year 
or from the corresponding Monte Carlo simulation samples. Values for the Mw and 
Fw estimators were determined separately for the data and simulation samples. This 
method has been applied on a cross-check analysis in the semi-leptonic channels and to 
the standard fully- hadronic analysis. The results from the fully-hadronic analysis are 
shown in Table 10. The semi-leptonic cross-check analysis applied the MLBZ procedure 
to the W mass determination separately in the electron, muon, and tau channels with 
uncertainties of around 8 MeV/c^ being obtained and the results being compatible with 
the systematic uncertainties quoted in this paper. The MLBZ method provides a useful 
cross-check of the size of the systematic uncertainty arising from fragmentation and other 
jet description errors reported in the previous section. Prom the values obtained from the 
MLBZ method we conclude that the systematic uncertainties have not been significantly 
underestimated. 

6.7 Elect roweak Radiative Corrections 

The measurements of the W mass and width described in this paper rely upon the 
accuracy of the event description provided by the simulation. Hence, the modelling 
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GeV 


AMw 
MeV/c^ 


MeV/c^ 


MLBZ 


qq'qq' Data - PYTHIA 


206.5 


-7.9 ± 4.9 


20.1 ± 10.5 



Table 10: Results obtained with the MLBZ method (see text). 



accuracy of the electroweak radiative corrections implemented in the event generator is 
a source of systematic uncertainty 

The radiative corrections for 4-fermion events are described in [13] and in section 4.2. 
For W+W^ (CC03) events, the signal used in this analysis, the corrections are based on 
YFSWW [17] and the effect of the theoretical uncertainties in it on the W mass measurement 
were initially studied in [34] at pure event generator level. 

In [35] this study has been performed in the context of the full DELPHI simulation 
and analysis procedure; furthermore the main uncertainties due to non-CC03 4-fermion 
background events have been studied. Radiative corrections uncertainties on non 4- 
fermion background events are included in the uncertainty estimated on the background. 

Several categories of uncertainty sources have been studied, which are considered here 
in turn. 

W+W" Production: Initial State Radiation (ISR) 

ISR plays a key role in the W mass analysis as it is one of the main sources of the 
bias on the fitted result with respect to the true value. This bias, which is removed by 
calibrating the fits with the simulation, is due to the energy-momentum conservation 
constraint used in the kinematical constrained fits. The ISR is computed in the YFS 
exponentiation approach, using a leading logarithm (LL) matrix element. 

The difference between the best result, obtained from implementing the O(a^) ISR 
matrix element, and the one provides an estimate of the effect of missing the 

matrix element for higher orders. The missing higher orders lead to the use of a wrong 
description for events with more than three hard photons or more than one photon with 
high pt. 

The difference between the best result and the 0{a) result includes the previous study, 
and can be used as an estimate of the upper limit of the effect of missing the non-leading 
logarithm (NLL) terms at (9(q;^); this effect of missing NLL terms is expected to be 
smaller than the effect from the LL terms given by this 0{a'^) to 0{a) difference. 

Also taking into account the study performed in [34], the ISR related uncertainty can 
be conservatively estimated at 1 MeV/ for the mass and 2 MeV/ on the width. 

W Decay: Final State Radiation (FSR) 

The FSR description and uncertainty is tightly linked to the final state considered. 
QED FSR from quarks is embedded in the parton shower describing the first phase of 
the hadronisation process. It is therefore essentially impossible to separate it from the 
rest of the hadronisation process, and the related uncertainty is considered as included 
in the jet and fragmentation related systematics. 

FSR from leptons is described by PHOTOS. The difference between the best result, 
based on the NLL treatment, and the LL one can give an estimate of the effect of the 
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missing part of the 0{a) FSR correction. While the result depends on the semi-leptonic 
channel, the difference is always less than 1 MeV/ c^. 

In [34] the effect of the missing higher orders beyond 0{a'^) has been found to be 
negligible at generator level. Simple perturbative QED considerations suggest that the 
size of the effect should not exceed the size of the effect from the missing part of the 
0{a) FSR correction; therefore conservatively the 1 MeV/ can be doubled to take into 
account both of these components of the uncertainty. 

Non-factorizable QED Interference: NF 0{a) Corrections 

Non-factorizable 0{a) QED interference between Ws is effectively implemented 
through the so-called Khoze-Chapovsky [36] (KC) ansatz. 

The effect of using the KC ansatz with respect to the Born calculation, where this 
interference is not described, can be considered as an upper limit of the missing part of 
the full 0{a) calculation and of the higher order terms. A dedicated study shows that 
the effect is less than 2 MeV/ for all the measurements. 

Ambiguities in Leading Pole Approximation (LPA) definition: Non Leading 
(NL) 0{a) Corrections 

Two sources of uncertainties are considered, following the study in [34]. The effect 
of missing higher orders can be, at least partly, evaluated by changing the electroweak 
scheme used in the 0{a) calculation. This essentially means changing the definition of 
the QED fine structure constant used in the 0{a) matrix element. The effect is very 
small, at the hmit of the fit sensitivity, both for the mass and the width. 

The second, more relevant, source of uncertainty connected to the LPA is in its possible 
definitions, i.e. the ambiguity present in the way of expanding the amplitude around the 
double resonant W pole. The standard YFSWW uses the so called LPA^ definition; a 
comparison with the LPA^ one can give an estimate of the effect from the intrinsic 
ambiguity in the LPA definition. A dedicated study has been performed evaluating the 
difference: 

AC(«)(LPA^-LPAb) = A(Best LPA^-no NL LPAA)-A(Best LPA^-no NL LPA^) 

in order to evaluate only the effect of the different scheme on the radiative corrections 
(and not at Born level). The size of the effect is less than 1 MeV/ for the mass and 
less than 4 MeV/c^ for the width. 

Radiative Corrections on 4-f Background Diagrams: Single W 

The Double Pole Approximation (DPA) is known to be valid within a few W widths 
of the double resonant pole. The DPA correction is applied only to the CC03 part of the 
matrix element (and partly to the interference, sec [13]); non-CC03 diagrams contribu- 
tions are not directly affected by the DPA uncertainty (except for possible effects in the 
interference term which is relevant for the electron channel). 

It is clear that this procedure still leaves the problem of the approximated radiative 
corrections treatment for the non-CC03 part of the matrix element (and the interference). 
The ISR studies previously discussed can reasonably cover the most relevant part of 
the electroweak radiative corrections uncertainties present also for the W"'"W~-like 4-/ 
background diagrams, e.g. the non-CC03 part. There is, however, a notable exception: 
the so called single W diagrams for the qq'eu final state. 

The bulk of single W events are rejected in the W mass and width analysis, since the 
electron in these events is lost in the beam pipe. But the CC03 - single W interference is 
sizeable, and it has a strong impact on the W mass result in the electron channel. The 
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situation is different in the W width analysis, where in eI7eqq' events reconstructed by the 
electron analysis the effects of non-CC03 diagrams and the CC03 - non-CC03 interference 
are opposite in sign and almost completely cancel. 

The situation is made even more complex by the cross-talk between channels, e.g. 
events belonging in reality to one channel but reconstructed as belonging to another one. 
This cross-talk is particularly relevant between semi-leptonic electron and tau decays, 
and this explains why the r channel analysis is also sensitive to this uncertainty source. 

The effect of this uncertainty has been studied in two ways. Firstly, since the uncer- 
tainty on the single W rate associated to radiative corrections is known in literature to 
be about 4%, the non-CC03 part of the matrix element, assumed to be dominated by the 
single W contribution, has been varied by 4% for qq'eu final states. Another possible 
source of uncertainty related to 4-/ background is estimated by partly applying the DPA 
correction to the interference term (see the discussion in [13]). The effect of this way 
of computing the corrections can be considered as another estimate of the uncertainty 
related to the 4-/ background presence. 

The maximal size of these effects is about 6 MeV/c^ (for the mass in qqeu and the 
width in qqrv). 

Total Uncertainty 

The results of all the studies presented are combined in a single uncertainty for each 
channel. Tables 11 and 12 present the estimates for the mass and width from the different 
sources of uncertainties discussed above. 



Mw Electroweak Correction Systematic Errors (McV/c^) 


Uncertainty Source 








qq'qq' 


ISR 


1.0 


1.0 


1.0 


1.0 


FSR 


0.5 


0.5 


1.0 




NF 0{a) 


1.0 


1.0 


1.0 


2.0 


NL 0{a) 


1.0 


1.0 


1.0 


1.0 


4-/ Background 


5.5 


0.5 


1.0 


0.5 


Total 


9 


4 


5 


4.5 



Table 11: Summary of the systematic uncertainties on the W mass due to electroweak 
corrections. The total is computed adding linearly the absolute values of all the contri- 
butions. 



Fw Electroweak Correction Systematic Errors (MeV/ c^) 


Uncertainty Source 


eI7eqq' 


/iz/^qq' 


rz7^qq' 


qq'qq' 


ISR 


2.0 


2.0 


2.0 


2.0 


FSR 


1.0 


1.0 


2.0 




NF 0{a) 


2.0 


2.0 


2.0 


2.0 


NL 0{a) 


4.0 


4.0 


4.0 


4.0 


4-/ Background 


2.0 


1.0 


6.0 


1.0 


Total 


11 


10 


16 


9 



Table 12: Summary of the systematic uncertainties on the W width due to electroweak 
corrections. The total is computed adding linearly the values of all the contributions. 
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Error [MeV] 


^/s Nominal [GeV] 


161 172 183 189 192 196 200 202 205 207 


25.4 27.4 20.3 21.6 21.6 23.2 23.7 23.7 36.9 41.7 



Table 13: Uncertainties on the Lep energies for the different centre-of-mass energy points. 

The total uncertainty per channel is conservatively computed summing linearly the 
values of the contributions. All the numbers have been rounded to 0.5 MeV/c^. 

Reference [13] also reports a comparison of YFSWW with the other completely indepen- 
dent Monte Carlo generator RacoonWW [37] which implements radiative corrections in the 
DPA. This study has not been directly used in the error estimation presented here due 
to the limitations in the treatment of non-coUinear radiation in RacoonWW. However, this 
study does provide additional confidence in the validity of the YFSWW calculation. 

As can be seen, the uncertainty on the W mass associated with the electroweak radia- 
tive corrections is found to be less than 10 MeV/ c^. 

6.8 LEP Collision Energy 

The average Lep collision energy is evaluated at 15 minute intervals of running or after 
significant changes in the beam energy. The measured centre-of-mass energy is imposed 
as a constraint in the kinematic fit, and hence the relative error on the collision energy 
translates to approximately the same fractional error on the W mass determination. The 
effect of the uncertainty on the W width determination is negligible. 

The beam energy is estimated using the Lep energy model, discussed in section 2 
based on 16 NMR probes in dipole magnets around the Lep ring calibrated with the 
RDP technique. The compatibility of three cross-check methods with this determination 
was used to determine a set of small energy offsets. The relative size of this offset was 
energy dependent, rising to a maximum of 1.6 x 10~^ at 207 GeV centre-of-mass energy. 

The Lep energy working group also assessed the uncertainties in the collision energies 
and supplied these in the form of a 10x10 correlation matrix. The uncertainties increase 
as the collision energy increases, due to the fact that higher energies are further from the 
RDP normalisation region. The errors are given in table 13. At 183 GcV centre- of- mass 
energy the uncertainty on the collision energy is 20.3 MeV. This rises to 23.7 MeV at 
202 GeV. For the energy points at values of 205 and 207 GeV, taken in the year 2000, 
there is an additional uncertainty due to the 'Bending Field Spreading' strategy, in which 
the corrector magnets were powered in a coherent manner to increase the overall dipole 
field and thus the Lep energy [7] . This leads to a larger error for the year 2000. For the 
energy points at 161 and 172 GeV, taken in the year 1996, there is also a small increase in 
the error, compared to 183 GeV, due to increased uncertainties in the NMR calibration 
for this year. 

The mean energy difference between the electron and positron beams is less than 
4 MeV at all energies and hence the effect on the W mass or width determination is 
negligible. The momentum spread of the electrons or positrons in a bunch gives rise 
to a variation in the centre-of-mass energy of the collisions and boost of the centre of 
mass frame with respect to the laboratory frame. The spreads in ccntrc-of-mass collision 
energies have been evaluated by the LEP energy working group [7] and range from 144 to 
265 MeV. The corresponding effects for the W mass and width analyses are negligible. 
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6.9 Aspect Ratio 

The aspect ratio is defined as the ratio of the length to the width of the detector. 
As all the sub-detectors of DELPHI are aligned with respect to the Vertex Detector, 
the knowledge of the aspect ratio is limited by the precision to which the position and 
dimensions of the Vertex Detector can be measured. The effect of a mismeasurement of 
the aspect ratio is to introduce a bias on the measurement of the polar angle, 6. As the 
W boson production polar angle is not isotropic but forward peaked, a mismeasurement 
of the aspect ratio would result in a small bias on the average opening angle of the W 
decay products, and hence induce a small bias on the reconstructed W mass. 

The correspondence of hits in the overlapping silicon modules is sensitive to a mis- 
alignment of the Vertex Detector. In fact the study of these overlaps constitutes an 
essential part of the procedure for the alignment of the Vertex Detector. From this study, 
discussed further in [9], it is concluded that a reasonable estimate of the aspect ratio 
uncertainty is 3 x 10~^. Such a bias would result in a shift in W mass below 1 MeV/ 
for the semi-leptonic channel, and of 2 MeV/ for the fuUy-hadronic one. The effect on 
the W width is negligible. 

6.10 Background Description 

The background events for the W-pair selection are from four-fermion or hadronic two 
fermion processes. 

The four-fermion background uncertainty is studied and described in the electroweak 
corrections uncertainties (section 6.7) and in the jet description studies (section 6.5) parts 
of this paper. 

The dominant source of background to W pair production, both in the semi-leptonic 
and in the fuUy-hadronic channel, is from Z — > qq{'~f) events. 

In the semi-leptonic channel the 2-fermion background is relatively small with the main 
uncertainty in its rate arising from the discrepancy between data and simulation in the 
rate of misidcntification of energetic photons (from radiative return to the Z peak events) 
as electrons. This misidcntification is mainly due to the electron-positron conversion 
of photons and the spurious associations of forward vertex detectors hits to an electro- 
magnetic cluster in the calorimeter. A data-simulation comparison shows that a 10% 
fluctuation of the background is possible without significantly degrading the agreement 
between the data and simulation. The theory uncertainty on the 2-fermion cross-section 
is generally small, in the worst case at the 2% level [38]. 

In the fuUy-hadronic channel the 2-fermion background is more important, and the 
major contribution to the uncertainty is from the four-jet final state production mecha- 
nism. The study performed in [39] has shown that the maximal difference in the estimated 
2-fermion background rate is 10% coming from changing from PYTHIA to HERWIG as the 
hadronisation model, with the ARIADNE model giving intermediate results. The effect on 
the W mass is 13 McV/c^ at = 189 GeV, and 4 MeV/c^ at ^/i = 206.5 GeV, while 
the effect on the W width is 40 MeV/ over the whole range of centre-of-mass energies. 

In summary, applying a variation of ±10% on the Z — > qq{j) event rate is used to 
provide an estimate of the systematic uncertainty on the background level for both the 
semi-leptonic and fuUy-hadronic channel mass and width measurements. This variation 
also covers any discrepancies seen in the data and simulation comparison plots shown in 
this paper. 

The importance of the background event mass distribution has also been investigated. 
In the semi-leptonic analyses the mass distribution taken from the simulation has been 
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replaced with a constant level and half of the variation in the result has been taken as a 
systematic. In the fuUy-hadronic channel this systematic was assessed by changing the 
generator used for the background between PYTHIA, HERWIG and ARIADNE. 

The background level and background shape uncertainties were added in quadrature 
and the resulting errors are reported in tables 14, 15 and 16 below. 

6.11 Bose-Einstein Correlations 

Correlations between final state hadronic particles are dominated by Bose-Einstein 
Correlations (BEC), a quantum mechanical effect which enhances the production of iden- 
tical bosons close in phase space. The net effect is that multiplets of identical bosons are 
produced with smaller energy-momentum differences than non-identical ones. 

BEC for particles produced from the same W boson affect the normal fragmentation 
and are therefore treated implicitly in the fragmentation uncertainties which are con- 
strained by the large amount of Z-data. BEC for pairs of particles coming from different 
Ws cannot be constrained or safely predicted by the information from single hadronically 
decaying vector bosons. 

A dedicated and model-independent measurement of the BEC effect was performed 
by the Delphi collaboration in [40] while other Lep experiments have made similar 
measurements [41]. Comparing these results with Monte Carlo models constitutes the 
only way to estimate potential systematic uncertainties from BEC. The LUBOEI model 
BE32 [42] was found to give the largest shift in the measured value of Mw for a given 
amount of BEC. Other models give smaller shifts and some models predict no appreciable 
BEC shifts at all. It was decided not to apply any corrections due to BEC and evaluate 
the systematic error as the largest predicted shift consistent with the Delphi data. The 
predicted shift plus one standard deviation of its error is used as the estimator of the 
systematic error. 

The Delphi result for BEC is a 2.4 standard deviation evidence for BEC between 
different Ws and a correlation strength, A, which can be compared to the BE32 prediction 
at the same effective correlation length scale: 

Adata/ABE32 = 0.55 ± 0.20(Stat.) ± 0.11(Syst.). (9) 

The predicted mass shift, BEC inside Ws only — BEC inside and between Ws, using 
BE32 (with model parameters PARJ(92) = 1.35 and PARJ(93) = 0.34) is 40±10 MeV/c^ 
for the standard mass analysis, 33±11 MeV/c^ for the cone jet mass reconstruction 
analysis and — 17±20 MeV/ for the W width analysis. The observed mass shift in BE32 
is linear in the observed correlation, Abe32- Applying the one standard deviation upper 
bound of the correlation parameter this translates into a systematic error of 31 MeV/ 
from BEC for the standard analysis and 26 MeV/ for the cone analysis. A systematic 
error of 20 MeV/ is apphed for the W width. The mass and width shifts were evaluated 
with the simulation model over the full range of centre-of-mass energies and no energy 
dependence was observed. The shifts reported are the average values. Conservatively, 
these errors are applied as symmetric uncertainties. 

The combined Delphi BEC measurements of the correlation strength and effective 
correlation length scale suggest that the between-W BEC occur with an effective corre- 
lation length scale which is larger that the one predicted by BE32. If this is the case, the 
number of pairs effectively affected by the BEC is reduced and also the effect per pair 
is diminished. Furthermore, the other Lep experiments have reported smaller values of 
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Adata/ Abe32 than that observed by DELPHI. Hence the systematic uncertainties apphed 
in this analysis are considered conservative. 

6.12 Colour Reconnection 

In the reaction e+e~ — * W^W~ (q'i^2)(q'3^4) the hadronisation models used for this 
analysis treat the colour singlets qiq2 and q^q^ coming from each W boson independently. 
However, interconnection effects between the products of the two W bosons may be 
expected since the lifetime of the W bosons (rw — fi/^w — 0.1 fm/c) is an order of 
magnitude smaller than the typical hadronisation times. 

The exchange of coloured gluons between partons from hadronic systems from different 
W bosons can induce the so-called colour reconnection (CR) effect in the development of 
the parton shower. This effect can in principle distort the properties of the final hadronic 
system and therefore affect the W mass measurement, if not properly accounted for in 
the simulation. 

At perturbative level the effects are expected to be small [43], and the impact on 
the reconstructed W mass has been evaluated to be at most 5 MeV/ c^. However, CR 
effects can be large at hadronisation level, due to the large numbers of soft gluons sharing 
the space-time region. These effects have been studied by introducing CR effects into 
hadronisation models and comparing with Delphi data and are reported in [30]. 

The most studied model, and the one used for the evaluation of the systematic uncer- 
tainty on the W mass and width measurement, is the Sjostrand-Khoze "Type 1" model 
(SK-I) [44]. This model of CR is based on the Lund string fragmentation phenomenol- 
ogy: the strings are considered as colour flux tubes with some volume, and reconnection 
occurs when these tubes overlap. The probability of reconnection in an event, P^eco) is 
parameterised by the value k, according to the volume of overlap between the two strings 

^overlap • 

Preco = 1 - 6-^^°^^^"^^. (10) 

The parameter k, determines the reconnection probability. By comparing the data with 
the model predictions evahiated at several k values it is possible to determine the value 
most consistent with the data and extract the corresponding reconnection probability. 

Another model has been developed by the same authors (SK-II') and also implemented 
in PYTHIA but is found to predict a smaller shift on the reconstructed W mass than SK-I 
for the same reconnection probability. 

Further CR models are available in the HERWIG and ARIADNE Monte Carlo programs. 
In ARIADNE, which implements an adapted version of the Gustafson-Hakkinen model [45] , 
the model used [46] allows for reconnections between partons originating in the same W 
boson, or from different W bosons if they have an energy smaller than the width of the W 
boson. The mass shift from CR is evaluated from the difference between the shift when 
the reconnections are made only in the same W boson and when the full reconnections 
are made. In the standard Delphi analysis, the shift was found to be 11 ± 11 MeV/ c^. 

In HERWIG the partons are reconnected, with a reconnection probability of 1/9, if the 
reconnection results in a smaller total cluster mass. The shift in the reconstructed W 
mass at 189 GeV centre-of-mass energy was found to be 29 ± 7 MeV/ c^, the same shift 
as obtained from a k value of 0.29 in the SK-I model. 

Delphi has performed two analyses to compare these simulation models with data 
which are described in detail in [30]. 

The first one is based on the measurement of the particle fiow between the jets in 
a four jets W+W~ event. On a subsample of strictly four-jet events two regions can 
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be defined, the region between jets from tlie same W (called inside-W regions) and the 
region between jets from different W bosons (called between- W regions). The ratio R of 
the particle fluxes in the inside-W and between-W regions (hmiting the analysis to the 
central part of these regions) is an observable sensitive to CR effects. The comparison of 
the flux measured in real data with the prediction of the SK-I model as a function of k 
allows the value to be determined which is most consistent with data, and its uncertainty. 

The second method used exploits the observation that in the direct reconstruction 
analysis of the W mass, different W mass estimators have different sensitivities to CR 
effects. As discussed in section 5.3.2 removing particles from the inter-jet regions reduces 
the sensitivity to CR effects and hence can be used to measure the CR effect. The 
correlation between the measurement of the mass shift (using the standard or cone jet 
reconstruction techniques) and the measurement of the mass from these techniques is 
only 11%. 

From the combination of these two analyses and in the framework of the SK-I model, 
the value of the k, parameter most compatible with the data is found to be [30] : 

K/ — 2.2 ^t-j^ 2 . 

The CR shift in the reconstructed W mass as a function of the SK-I k parameter is 
provided as figure 13, the results of the standard and cone jet reconstruction techniques 
are indicated. Figure 14 shows the CR shift for the W width reconstruction analysis. 

The systematic uncertainty on the W mass and width is calculated using the one 
standard deviation upper bound of k, of 4.7. As reported above, this systematic error 
is considerably larger than that which would be evaluated from the ARIADNE or HERWIG 
CR models. Furthermore, this value of k is larger than that reported by the other Lep 
experiments [31]. The CR W mass shift is dependent on the centre-of-mass energy in the 
SK-I model as shown in figures 13 and 14. However, we prefer not to rely on the centre- 
of-mass energy evolution of the SK-I CR shift (leading to a change in relative weights 
when averaging the results from different centre-of-mass energies) and instead choose to 
quote the systematic errors at 200 GeV (close to the average centre-of-mass energy of the 
data). In light of the significant range of CR effect estimates no correction is made to the 
W mass or width results and for simplicity a symmetric systematic uncertainty is applied. 
The corresponding systematics uncertainties on the W mass are 212 MeV/ (standard), 
116 MeV/ (cone jet reconstruction) and 247 MeV/ for the W width analysis. 

7 Results 

The results of the analyses and the final combinations of these results are presented in 
this section. The results are obtained at a range of nominal centre-of-mass energies and 
in the four event selection channels. Combined results are obtained from an average of 
these results and also an average with the previously published Delphi data [1,2] that 
have not been reanalysed in this paper. 

Subdividing the results by data-taking years and nominal centre-of-mass energies en- 
ables a proper treatment of the correlated systematic uncertainty from the Lep collision 
energy and other dependences on the centre-of-mass energy or data-taking period. A 
detailed breakdown of the sources of systematic uncertainty, as shown in tables 14,15 and 
16, is provided for each result and the correlations specified. 

The combination is performed and the evaluation of the components of the total error 
assessed using the Best Linear Unbiased Estimate (BLUE) technique [47]. 
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Figure 13: W mass shift caused by the colour reconnection effect as described in the 
SK-I model plotted as a function of the model parameter k which controls the fraction 
of reconnected events. The upper plot is for the standard W mass analysis and the lower 
plot when the cone jet reconstruction technique is applied. 
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Figure 14: W width shift caused by the colour reconnection effect as described in the 
SK-I model plotted as a function of the model parameter k which controls the fraction 
of reconnected events. 
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Mw Systematic Errors (MeV/c^) at 189 GeV 



oouiceb oi oybieiiiaiic njiioi 






Ti^T-qq' 




Statistical Error on Calibration 


12 


10 


15 


4 


Lepton Corrections 


24 


16 






Jet Corrections 


18 


15 


19 


24 


Fragmentation 


10 


10 


13 


12 


Electroweak Corrections 


9 


4 


5 


5 


Background 


5 


1 


12 


17 


LEP Energy 


9 


9 


9 


9 


Bose-Einstein Correlations 

Colour Reconiiection 




31/26 
212/116 



Table 14: Contributions to the systematic error on the W mass measurement for data 
taken at a nominal centre-of-mass energy of 189 GeV. Where two uncertainties are 
reported in the qq'qq' analysis column the first corresponds to the standard analysis and 
the second to the cone jet reconstruction analysis. 

7.1 W Mass 

The W mass is extracted separately in the analyses designed to select the el^eqq', 
//I/^qq' and ri^^qq' decay channels. The values obtained are given in table 17 for the 

analysed centre-of-mass collision energies. The semi-leptonic channel analysis results are 
combined into a single fl/^qq' value for each year of data taking. When performing these 
combinations the following sources of systematic uncertainty are taken as fully-correlated 
between lepton channels and between years: electroweak corrections, fragmentation, jet 
corrections, lepton corrections, background. The Lep energy mcasTircmcnt correlations 
are taken from the matrix supplied in [7]. The simulation calibration statistics are taken 
as uncorrelated. 

The W mass is also obtained from the qq'qq' channel using both the standard and 
cone jet reconstruction technique. The results obtained from these analyses are given in 
table 18. 

In addition to the analyses presented in this paper, measurements of the W mass have 
also been made using the data collected in 1996. 

7.1.1 W Mass from the W+W~ Cross-section 

The Delphi collaboration has measured the total CC03 W+W" cross-section, as a 
function of centre-of-mass energy, using the full data sample collected by the collabora- 
tion during Lep2 operations [39]. Assuming the validity of the cross-section dependence 
predicted by the Standard Model these measurements can be translated into a measure- 
ment of the W mass. Only the cross-section measurements close to the W+W~ threshold 
have significant sensitivity to the W mass. 

The Standard Model cross-section dependence on the W mass is obtained from the 
WPHACT and YFSWW generator setup, as discussed in section 4.2, and cross-checked with 
the improved Born approximation calculation. The theoretical error on the total W+W~ 
cross-section near threshold was estimated as 2% decreasing with increasing collision 
energy to 0.5% in the DPA- valid region [48], the corresponding error on the W mass is 
marked below as Theor. The sources of experimental systematic error have not been 
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Mw Systematic Errors (MeV/c^) at 205 GeV 
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Table 15: Contributions to the systematic error on the W mass measurement for data 
taken at a nominal centre-of-mass energy of 205 GeV. Where two uncertainties are 
reported in the qq'qq' analysis column the first corresponds to the standard analysis and 
the second to the cone jet reconstruction analysis. 
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Table 16: Contributions to the systematic error on the W width measurement for data 
taken at a nominal centre-of-mass energy of 205 GeV. 

reevaluated and are as reported in [1], apart from use of the revised collision energy 
uncertainty. 

From a fit of the measured cross-sections at centre-of-mass energies of 161.31, 172.14 
and 182.65 GeV the mass has been determined to be 

Mw = 80.448 ± 0.434(Stat.) ± 0.090(Syst.) ± 0.043(Theor.) ± 0.013(LEP) GeV/c^ 

7.1.2 W Mciss from Direct Reconstruction at = 172 GeV 

For completeness, we also report here on the relatively small data sample (10 pb"-*^) 
recorded in 1996 at ^/s = 172 GeV. This sample was analysed and W mass results 
published using the eI7eqq', A^i^/^qq' and qq'qq' decay channels in [2]. The qq'qq' analysis 
was performed using a standard analysis rather than a cone jet reconstruction based 
analysis. 

This data sample has not been reprocessed, nor have W width results been produced 
with this sample. The estimates of systematic uncertainties are retained from the original 
paper except for the uncertainties arising from colour reconnection and Bose-Einstein 
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Correlations in the qq'qq' channel, where the errors reported above for the standard 
analysis are used, and the use of the final Lep collision energy uncertainty. The revised 
values are 

Mw = 80.51 ± 0.57(Stat.) ±0.05(Syst.) ±0.01(LEP) GeV/c^ 

for the combined semi-leptonic channels, and 

Mw = 79.90 ± 0.59(Stat.) ± 0.05(Syst.) ± 0.21(FSI.) ± O.Ol(LEP) GeV/c^ 

for the fuUy-hadronic decay channel. These values have been included in tables 17 and 
18. 

7.1.3 Combined Results 

The combinations of the results arc performed, assuming that the following compo- 
nents of the error are fully-correlated between years (and energy points) and between the 
fuUy-hadronic and semi-leptonic channels: electroweak corrections, fragmentation and 
jet correction. The lepton-related detector systematic in the semi-leptonic channel is 
also assumed to be fully correlated between years. The colour reconnection and Bose- 
Einstein effect in the fully-hadronic channel is assumed to be fully correlated between 
years. The error arising from calibration statistics is uncorrelated between years in the 
semi-leptonic analysis, as it was determined from independent Monte Carlo simulation 
samples, but this error is correlated in the fully-hadronic channel as the values were ob- 
tained from an overall fit to the samples at all centre-of-mass energies. This error source 
is uncorrelated in the combination of the semi-leptonic and fully-hadronic channel. The 
background-related systematic is assumed to be fully correlated between years in both 
the fully-hadronic and semi-leptonic analyses but uncorrelated between the two channels. 
The LEP centre-of-mass energy uncertainty is, of course, fully correlated between the 
semi-leptonic and fully-hadronic decay channels but is only partially correlated between 
years. The inter- year correlations were assessed by the LEP energy working group [7] and 
this correlation matrix was applied when performing the combinations reported here. 

The results from the semi-leptonic W mass analyses in each year of data taking (1996- 
2000) have been combined. The result for the analysis aimed at selecting events in the 
ei/gqq' decay channel is: 

Mw = 80.388 ± 0.133(Stat.) ± 0.036(Syst.) ± O.OIO(LEP) GeV/c^ 

the combination has a probability of 25%. 

The result for the analysis aimed at selecting events in the //i/^qq' decay channel is: 

Mw = 80.294 ± 0.098(Stat.) ± 0.028(Syst.) ± O.OIO(LEP) GeV/c^ 

the combination has a probability of 96%. 

The rz7^qq' selection includes significant cross-talk from events in other decay channels 
(see table 2) and a result from the 1996 data is not available. The result for the analysis 
aimed at selecting events in the rl/^qq' decay channel (in the years 1997-2000) is: 

Mw = 80.387 ± 0.144(Stat.) ± 0.033(Syst.) ± O.OIO(LEP) GeV/c^ 

the combination has a probability of 56%. 

The result for the combined semi-leptonic W mass analyses is: 

Mw = 80.339 ± 0.069(Stat.) ± 0.029(Syst.) ± 0.009(LEP) GeV/c^ 
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the combination has a probabihty of 16%. 

Similarly, the results on the W mass extracted from the fuUy-hadronic event analysis 
have also been combined. The value from 1996 uses the standard reconstruction tech- 
nique; the results of the cone-jet reconstruction technique are used for the other data 
taking years (1997-2000). The combined result is: 

Mw = 80.311 ± 0.059(Stat.) ± 0.032(Syst.) ± 0.119(FSI) ± O.OIO(LEP) GeV/c^ 

the combination also has a probabihty of 16%. 

The mass difference between the W boson mass measurements obtained from the 
fuUy-hadronic and semi-leptonic channels AMw(qq'qq' — ^i^^qq'); has been determined. 
A significant non-zero value for AMw could indicate that Bose- Einstein or colour recon- 
nection effects are biasing the value of Mw determined from qq'qq' events. Since AMw 
is primarily of interest as a cross-check of the possible effects of final state interactions, 
the errors from CR and BEG are set to zero in its determination and the results of the 
standard reconstruction technique, rather than the FSl effect-reducing cone-jet recon- 
struction technique, are used for the qq'qq' analysis. The result provides no evidence for 
FSI effects: 

AMw(qq'qq' - ^I^^qq') = 0.024 ± 0.090 GeV/c^, 

the combination has a x^ probability of 20%. 

The final DELPHI result for the W mass for the full Lep2 data sample is obtained 
by combining the values obtained from the direct reconstruction method in the iViqc^ 
analysis and cone jet reconstruction technique qq'qq' analysis in each data taking year. 
The value obtained from the threshold cross-section is also included in this average. The 
combined result is: 

Mw = 80.336 ± 0.055(Stat.) ± 0.028(Syst.) ± 0.025(FSI) ± 0.009(LEP) GeV/c^ 

the combination has a x^ probability of 15%. 

Although the statistical error in the iP^qq' and qq'qq' channels is similar, owing to 
the large systematic error attributed to final state cross-talk effects the weight of the 
fuUy-hadronic channel results in this average is 21%. The weight of the threshold cross- 
section measurement of the W mass is only 2% due to the small data sample collected at 
161 GeV centre-of-mass energy. The full error breakdown of the averages is provided in 
table 19. 

The Delphi measurement of the colour reconnection effect is reported in [30]. This 
measurement places relatively loose constraints on the size of the W mass uncertainty from 
GR effects, and thus leads to the small impact of the fuUy-hadronic mass in the Delphi 
average. For comparison the value of the combined Delphi W mass as a function of 
the GR uncertainty is shown in table 20. All other errors, including that arising from 
Bose-Einstein correlations, have been kept constant in these results. 

7.2 W Width 

The W width has been measured from the semi-leptonic and the fuUy-hadronic decay 
channel events. As the analysis is less sensitive to the W width than the W mass, the 
width is extracted by performing a combined fit of the three semi-leptonic channels rather 
than from each channel individually. The results are given in table 21. The correlations 
assumed for the combinations are identical to those reported above for the W mass. 
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The results from the semi-leptonic W width analyses in each year of data taking 
(1997-2000) have been combined, the result obtained is: 

Tw = 2.452 ± 0.184(Stat.) ±0.073(Syst.) GeV/c^ 

the combination has a probability of 9%. 

Similarly, the results on the W width extracted from the fuUy-hadronic event analysis 
have also been combined, the result obtained is: 

Tw = 2.237 ± 0.137(Stat.) ± 0.139(Syst.) ±0.248(FSI) GeV/c^ 

the combination has a probability of 62%. 

The final Delphi result for the W width for the full Lep2 data sample is obtained 
by combining the values obtained from the direct reconstruction method in the £l7^qq' 
analysis and qq'qq' analysis in each data taking year. The combined result is: 

Tw = 2.404 ± 0.140(Stat.) ±0.077(Syst.) ±0.065(FSI) GeV/c^ 

the combination has a probability of 27%. 

Although the statistical error in the fz/^qq' and qq'qq' channels is similar, owing to 
the large systematic error attributed to final state cross-talk effects the weight of the 
fuUy-hadronic channel results in this average is 26%. The full error breakdown of the 
averages is provided in table 22. 

8 Conclusions 

The mass and width of the W boson have been measured using the reconstructed 

masses in e"'"e~^W"'"W~ events decaying to qq'qq' and fz7^qq' states. The W Mass was 
also extracted from the dependence of the W"^W~ cross-section close to the production 
threshold. The full Lep2 data sample of 660 pb^^ collected by the Delphi experiment 
at centre-of-mass energies from 161 to 209 GeV has been used. The final results are: 

Mw = 80.336 ± 0.055(Stat.) ± 0.028(Syst.) ± 0.025(FSI) ± 0.009(LEP) GeV/c^ 

Tw = 2.404 ± 0.140(Stat.) ±0.077(Syst.) ±0.065(FSI) GeV/cl 
These results supersede the previously published Delphi results [1-4]. 
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Table 17: Measured W mass (in GeV/c^) from the semi-leptonic decay channel analyses 
with the nominal centre-of-mass energies (in GeV) of each data sample indicated. The 
values marked iV^qc^ are the combined values of the three semi-leptonic channel analyses. 
The values obtained from the data recorded in 1996 and analysed in [2] are also included. 
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Year 
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80.318 ± 
80.171 ± 



0.092(Stat. 
0.104(Stat. 



± 0.032(Syst. 
± 0.032(Syst. 



± 
± 



0.214(FSI 
0.119(FSI 



± 0.015(LEP) 
± 0.015(LEP) 



Table 18: Measured W mass (in GeV/c^) from the fuUy-hadronic decay channel analysis 
with the nominal centre-of-mass energies (in GeV) of each data sample indicated. Results 
are provided for both the standard (std) and cone jet reconstruction techniques applied. 
The value obtained from the data recorded in 1996 and analysed in [2] is also included. 







qq'qq' 


All 


Value 


80.339 


80.311 


80.336 


Statistical Error 


.069 


.059 


.055 


Statistical Error on Calibration 


.003 


.004 


.002 


Lepton Corrections 


.015 




.012 


Jet Corrections 


.020 


.026 


.021 


Fragmentation 


.011 


.012 


.011 


Background 


.007 


.013 


.006 


Threshold Systematics 






.002 


Electrowcak Corrections 


.006 


.005 


.006 


LEP Energy 


.009 


.010 


.009 


Bose-Einstein Correlations 




.026 


.005 


Colour Reconnection 




.116 


.024 



Table 19: The final results (in GeV/ c^) of the W mass analyses and the breakdown of the 
uncertainty into its component categories. The ^^qq' and qq'qq' results use the values 
obtained in these analysis channels from the direct reconstruction method. The column 
marked 'All' uses the full direct reconstruction analyses and the threshold cross-section 
measurement. The qq'qq' results are taken from the cone jet reconstruction analysis, for 
all data except 1996 where the standard analysis was used. 



59 



CR MeV/c^ 



'^SK-I 



Mw GeV/c^ 





20 
40 
60 

80 
100 



0.00 
0.40 
0.89 
1.51 
2.30 
3.36 



80.326 ± 0.045(Stat.) 
80.326 ± 0.045(Stat.) 
80.328 ± 0.046(Stat.) 

80.330 d= 0.048(Stat.) 
80.333 ± 0.051(Stat.) 
80.335 ± 0.054(Stat.) 



± 0.028(Syst. 
± 0.028(Syst. 
± 0.028(Syst. 
± 0.028(Syst. 
± 0.028(Syst. 
± 0.028(Syst. 



) ± 0.013(FSI) 
) ± 0.016(FSI) 
) ± 0.021(FSI) 
) ± 0.024(FSI) 
) ± 0.026(FSI) 
) ± 0.026(FSI) 



± O.OIO(LEP) 
± O.OIO(LEP) 
± O.OIO(LEP) 
± O.OIO(LEP) 
± O.OIO(LEP) 
± 0.009(LEP) 



Table 20: The combined DELPHI W Mass value as a function of the uncertainty ascribed 
to colour reconnection effects in the fuUy-hadronic decay channel. The values of the 
'^SK-I parameter that give rise to this shift in the qq'qq' W mass at a centre-of-mass 
energy of 200 GeV are also given. 



Year 


Energy 


Channel 


Fw GeV/c^ 


1997 


183 


£i7^qq' 


2.495 ± 0.590(Stat.) 


± 0.069(Syst.) 








qq'qq' 


2.572 ± 0.460(Stat.) 


± 0.092(Syst.) 


± 0.248(FSI) 


1998 


189 


£z7^qq' 


3.056 ± 0.401(Stat.) 


± 0.071(Syst.) 








qq'qq' 


2.337 ± 0.260(Stat.) 


± 0.114(Syst.) 


± 0.248(FSI) 


1999 


192 


£z7^qq' 


2.342 ± 0.953(Stat.) 


± 0.071(Syst.) 








qq'qq' 


2.390 ± 0.756(Stat.) 


± 0.126(Syst.) 


± 0.248(FSI) 




196 


£z7^qq' 


1.805 ± 0.440(Stat.) 


± 0.072(Syst.) 








qq'qq' 


2.545 ± 0.508(Stat.) 


± 0.142(Syst.) 


± 0.248(FSI) 




200 


£z7^qq' 


2.153 ± 0.477(Stat.) 


± 0.073(Syst.) 








qq'qq' 


2.210 ± 0.376(Stat.) 


± 0.157(Syst.) 


± 0.248(FS1) 




202 


£z7£qq' 


1.707 ± 0.649(Stat.) 


± 0.076(Syst.) 








qq'qq' 


1.797 ± 0.488(Stat.) 


± 0.165(Syst.) 


± 0.248(FSI) 




192-202 




1.950 ± 0.277(Stat.) 


± 0.072(Syst.) 








qq'qq' 


2.210 ± 0.243(Stat.) 


± 0.152(Syst.) 


± 0.248(FSI) 


2000 


206 


£l7^qq' 


2.814 ± 0.364(Stat.) 


± 0.083(Syst.) 








qq'qq' 


1.979 ± 0.225(Stat.) 


± 0.183(Syst.) 


± 0.248(FSI) 



Table 21: Measured W widths (in GeV/c^) from the semi-lcptonic decay and ful- 
ly-hadronic decay channel analyses with the nominal centre-of-mass energies (in GeV) of 
each data sample indicated. 



60 







qq'qq' 


All 


Value 


2.452 


2.237 


2.404 


Statistical Error 


.184 


.137 


.140 


Statistical Error on Calibration 


.006 


.009 


.005 


Lepton Corrections 


.041 




.030 


Jet Corrections 


.036 


.129 


.059 


Fragmentation 


.029 


.008 


.024 


Electroweak Corrections 


.011 


.009 


.010 


Background 


.037 


.051 


.031 


Bose-Einstein Correlations 




.020 


.005 


Colour Reconnection 




.247 


.065 



Table 22: The final results (in GeV/c^) of the W width analyses and the breakdown 
of the uncertainty into its component categories. The -^I^^qq' and qq'qq' results use the 
values obtained in these analysis channels from the direct reconstruction method. The 
column marked 'All' provides the result from combining the measurements made in both 
channels. 



